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FREE-STREAM BOUNDARIES OF TURBULENT FLOWS ^ 


By Stanley Corrsin and Alan L. Kistler 


SUMMARY 

An experimental and iheoreiical study has been made of the 
instantaneously sharj) and irregular front which is always 
found to separate turbulent fluid from contiguous ^^nonturbu- 
l end ^ fluid at a free-stream boundary. This distinct demarca- 
tion is known to give an intermittent character to hot-wire 
signals in the boundary zone. 

The overall behavior of the front is described statistically in 
terms of its wrinkle-amplitude growth and its lateral propaga tion 
relative to the fluid as functions of downstrearn coordinate. 

It is proposed and justified that the front actually consists 
of a very thin fluid layer in which direct viscous forces j)lay the 
central role of transmitting mean and fluctuating vorticity to 
previously nonturbulent fluid. Outside this ^‘laminar super- 
lay eT^ there is presumably a field of ir rotational velocity 
fluctuations {the ^^nonturbulenV^ flow) with constant mean 
velocity. As outlined in the following j^aragraphs, theoretical 
analysis based on this general physical picture gives results on 
front behavior which are in plausible agreement with experi- 
mental results for three turbulent shear flows: rough-wall 
boundary layer ^ plane wake, and round jet. 

It is shown that the rate of increase of wrinkle amplitude of 
the front can be roughly exjdained as a Lagrangian diffusion 
process, using the statistical properties of the turbulence in the 
fully turbulent zone. 

The transversal propagation velocity of the turbulence front 
is predicted by the behavior of a physicomathematical model of 
the laminar superlayer. The model is a generalized Stokes- 
Rayleigh infinite wall, oscillating in its own plane, translating 
to give constant mean vorticity at the boundary, local 

vorticity production and uniform suction velocity. 

Finally, various statistical properties of the turbulence front 
location as a stationary random variable {for fixed downstream, 
position) have been either directly measured or indirectly inferred 
from known theorems on Gaussian stochastic processes; it is 
found that for boundary layer, wake, and jet the front location 
is very nearly Gaussian. Speciflcally, it is possible, therefore, 
to estimate the autocorrelation function of the front position, 

INTRODUCTION 

Until the last few years, basic experimental and (especially) 
theoretical attacks upon the problems of turbulent flow have 


centered on fully turbulent fields, both isotropic and shearing. 
The experimental researches have been concerned with the 
measurement of significant statistical quantities with the 
hope that these will give some insight into the mechanism of 
fully developed turbulence and might even suggest a profit- 
able theoretical approach. 

In reality, however, every turbulent flow is bounded l)y 
fluid not in a turbulent state. If the boundary spacings can 
be made veiy large compared with the characteristic correla- 
tion lengths of the turbulence, for example, integral scale and 
dissipative scale (microscale), then an “infinite field” approxi- 
mation can be used. This has been possible in research on 
the decaying turbulence behind regular grids, a reasonably 
good likeness of Taylor’s ideal concept of isotropic turbulence 
(ref. 1). 

It now seems probable that the classic turbulent shear 
Hows, boundary layer, wake, jet, channel, and so forth have 
transversal integral scales not very small compared with 
their characteristic widths. This has been shown experi- 
mentall}^ for the round jet (ref. 2), the plane half jet (ref. 3), 
the boundary layer (ref. 4), and the channel (ref. 5). This 
implies that the general behavior of these shear flows cannot 
be fully inferred on a (still unsolved) homogeneous shear flow^ 
basis but must involve the boundary phenomena. 

Turbulent shear flow^ boundaries can be classified in various 
ways. A conventional one is the division into (a) solid and 
(b) free (or free stream) boundaries, depending upon the 
presence or absence of a solid wall and excluding possible 
symmetry planes from consideration as boundaries. 

A further subdivision can be made in each class according 
to whether the outside flow' or wall is traveling faster or 
slower than the turbulent fluid just inside the boundary, 
but this distinction is probably only a quantitative one 
(because of the nonlinearity of the sy^stem), not affecting the 
nature of the boundary phenomena; a comparison of wake 
and jet boundaries would illustrate this remark. One can 
also visualize a boundary state in which this mean velocity 
difference is zero, that is, the case of uniform velocity field 
including both turbulent and outside flow*. 

This investigation is concerned solely with the free 
boundary condition. In practice, this case generally involves 
a mean shear stress in the fully turbulent region, reducing 


1 Supersedes NACA TN 313.3, “The Free-Stream Boundaries of Turbulent Flows” by Stanley Corrsin and Alan L. Kistler, 1954. 
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to zoro inonotonicall}^ toward the iionturhulent ^ free-stream 
flow. 

The outstanding observal)le characteristic of free bound- 
aries is tlie relativel^^ sharp instantaneous demarcation 
surface l)etween turbulent and nonturbulent fluid. This 
shows iij) very clearly, for example, in short-duration shadow- 
gra])hs of the turbulent wakes behind high-speed projectiles 
(fig. 1). The sharpness of the irregular boundary illustrated 
persists as far downstream as pictures have been taken, 
about several hundred wake diameters. 



L-81230 

Figure 1.— Turbulent wake of bullet. (Courtesy of Ballistic Research 
Laboratories, Aberdeen Proving Ground.) 

In a mixed flow zone of this type, a probe stationary rela- 
tive to the disturbance (e. g., the wall in a turbulent boundary 
layer) will be swept over by successive sections of turbulent 
and nonturbulent fluid. ith a hot-wire anemometer this 
Auelds an intermittent signal of the type which led to the 
discovery of this characteristic of free turbulent boundaries 
(ref. 2). The relative time spent by the probe in turbulent 
fluid was first measured b}" Townsend (ref. 0) and called 
the intermittency factor 7. 

Most steady-state shear zones spread with increasing 
downstream distance. Therefore, there cannot be even 
rough overall flow similarity unless the average lateral 
position and the wrinkle amiilitude of the sharj) boundary 
both increase at roughly the same rate as does the momentum 
width of the shear flow. Since it is well known that most 
‘‘simple’’ turbulent shear flows exhibit a rough overall 
similarity, it can immediately be anticipated that this 
turbulence front must (a) jiropagate relative to the local 
fluid in the same sense liiat a flame front propagates through 
a combustible mixture and (b) increase its geometrical 
amplitude with increasing downstream coordinate. 

The explanations of these necessary jiroperties of the 
turbulence front are two of the explicit ])urposes of this 
investigation. The two jiroiierties are to be measured and 
to be analytically related to jdiysical properties of the turbu- 
lence in the full}' turbulent zone. 

For any a:-station, the intermittency factor 7(2/) is just 1 
minus the distribution function of }^t), the instantaneous 
location of the sharp front between turbulent and nonturbu- 


2 The term “laminar” is reserved for a nonturbulent flow in shear, that is, where viscous 
forces are important. This is in contrast with the terminology introduced in reference 2, 
where laminar was used to indicate any nonturbulent flow. Of course, in practice, a “nou- 
turbulent" flow may be one whose turbulence level is much lower than that of the con- 
tiguous turbulent flow. 


lent fluid. For a fixed value of j:, Y(f) is a stationary 
random variable, and 

y(y) =i)i*ob[?/^ V{t) ^00] (1) 


Since y(y) is differentiable (in fact, nondiffereiitiable func- 
tions cannot be experimentally so identified), dyldij is the 
])robability density of ¥{(). 

A priori the fact that the free turbulence boundary 
(vorticity fluctuation boundary according to the physical 
])icture ])roposed here) remains sharp can be attributed to 
the continuous irregular stretching of the local vorticity 
gradient in the boundary, that is, to the fact that the vorticity 
])ro])agation process is nonlinear; for a given stretching rate, 
the production of new vorticity is jiroportional to the amount 
already present. This must be balanced on the average 
by the viscous diffusion of the vorticity gradient at the front. 

It is obvious that the random vorticity field ordinarily 
called turbulence can propagate only by direct contact, as 
opjiosed to action at a distance, because rotation can be 
transmitted (o irrotatiorial flow only through direct viscous 
shearing action. This insures that under ordinary circum- 
stances the turbulence front will always be a continuous 
surface; there will be no islands of turbulence out in the free 
stream disconnected from the main body of turbulent 
fluid. 

The anal}tical estimates will include a hyjiothetical case 
in which tlu‘ tuiFulent part of the flow field is also witliout 
shear. This is perhajis the simplest conceivable case under 
which turbulence ])ropagates into nonturbulent fluid — 
])rovided that one can neglect the necessary monotonic time 
decreases in turlndeiit energy ])er unit mass. Under these 
conditions it is projiosed that the distinction lietween turbu- 
lent and nonturbulent zones is the ])resence or absence, 
respectively, of random vorticity fluctuations. 

A more complex case is the one ordinarily encountered in 
])ractice, as described before: a shearing turbulence en- 
croaching on a nonshearing (irrotational) nonturbulent fluid. 
In (his case, the average propagation velocity of the turbu- 
lence front should also depend upon the mean shear stress 
in the tuiFulent fluid near the front. 

A somewhat different situtation, not included fully in the 
above classes, occurs in the transitional spreading of a turbu- 
lent shear region into a shearing laminar region, when the 
principal shear planes of laminar and turbulent flows are 
jiarallel to each othei- but perj)eiidicular to the mean propa- 
gation front. 

Such a phenomenon was first studied experimentally by 
Charters (ref. 7 ), who called it “transition by transverse 
contamination.’^ Emmons (ref. 8) has given good expei’i- 
mental evidence that transition from laminar to turbulent 
flow may often occur in this way, usually from ircegularly 
generated ‘•ignition” spots in the moving fluid, and a pre- 
liminary analytical discussion of the turbulence spread under 
these conditions has been given by Mitchner (ref. 9 ^ How- 
ever, it apjrears that Mitchner has omitted from his non- 
turbulent ri'gion the very shear which distinguishes the tran- 
sition problem. It is not intended that this important case 
l)e included in the present report. Although some of the 
same phenomena may occur as in the simpler nonshear 
boundary, it is possible that the dominant turbulence propa- 
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gatioii medumism is different. In particular, it may be 
that a destabilization of tbc already rotational flow occurs in 
addition to a transmission of random vorticity by direct 
viscous action at the turbulent-laminar interface. 

AVlien the present work was begun, it was hoped that the 
problem of propagation of turbulence into a nonturbulent 
flow could be studied at the boundary between a grid- 
generated isotropic turbulence and a nonturbulent flow 
moving at the same uniform speed. This would eliminate 
the shear stress entirely, although involving a relatively ra])id 
turbulence-level change due to viscous dissipation with no 
turbulent energy production. 

The principal generating arrangement tried was a half grid 
consisting of a conventional 1 -inch-mesh, )i-inch-dowel grid 
covering half the tunnel cross section, with a fine mesh 
screen of virtually identical static-pressure drop covering 
the other half. Unfortunately, anomalous boundary be- 
havior, arising from complexities in the flow around the joint 
between grid and screen, could not be eliminated with a 
reasonable amount of eflbrt. I'herefore, the turbulence 
propagation has been studied in situ, chiefly at the outer 
edge of a low-speed turbulent boundary layer, with a few 
measurements in a round jet for an additional check of some 
particular phenomena. For completeness, some of Town- 
semrs plane-wake data (ref. 10) have also been analyzed in 
the light of this investigation. 

The general purpose of this investigation has been to 
measure statistical projierties of the propagating turbulence 
front to permit qualitative or even rough quantitative theo- 
retical explanation of the phenomenon. 

''fhe work has been carried out at the Department of 
Aeronautics of The Johns Hopkins University with the 
financial assistance and sponsorship of the Xational Advisory 
Committee for Aeronautics. The authors would like to 
acknowledge the assistance of Miss Vivian O’Brien, Mr. 
Aristoteles Scoledes, Mr. Donald Johnson, and Miss M. Ann 
Emmart as well as the critical advice of Dr. Francis II. 
Clauser and Dr. Mark V. Morkovin. 


L 

U 

Az, 

luU 

M 

m 

Nr 

Ny, No 


n 

P 

PU P2 


0 

g 

9 . 

R 

Pl 

Pu. 

Pe, l/A, 2/A 
Px 
r 
r 

o 


S 


s=S-S 
T, T, ■ 


SYMBOLS 


A 

a 


Cly 

B 

b 

r, I) 

C/ 


a 

Fc, Fa, Fa^, Fa 


r 


characteristic ordinate in sketch (d) 
random variable re]>resenting some How 
]n-operty 

same property, taken in turbulent flow only 
constant 

random on-off signal, taken between zero 
and 1 

random variables 
skin friction coefficient, 



diameter of rod used to produce plane wake 

power spectra 

scalar function 

height of wall roughness 

random variable 

parameter in model of laminar superlayer, 
equivalent to vortex stretching rate 
wave number 


t 

u, 1^ ir 

Vo 

v„ 

i\ 

u,v,w 

Ui = ([ 

U* " 






Xi = r 


transversal Euler ian scale 

Lagrangian time scale 

Lagrangian length scale, v'Lt 

average pulse lengths of intermittent signal 

empirical constant 

exponent of boundary-la^^er j)ower-law ve- 
locity profile 

average frequency of occurrence of any 
particular value of random variable I{t) 
average frequency of occurrence of Y and 
zero, respectively, in front-location vari- 
al)le Y{t) 
cyclic frequency 
total static pressure 

probability densities of turbulent and po- 
tential segment lengths, respectively, in 
intermittent signal 
total velocity vector 
velocit^^ (luctuation vector, 
dynamic pressure in free stream of wind 
tunnel 

instantaneous radial location of turbulence 
front in round jet 
Lagrangian correlation function 
shear correlation coefficient, tTy/u'y' 
Reynolds numbers of laminar superlayer 
turbulence Reimolds number, 
coordinate vector 
radial coordinate in round jet 
jet orifice radius 

radial position at which U Umax at a sec- 
tion of the jet 

total shear force vector (per unit area) at 
turbulent side of superlayer, lying in 
plane of supei-layer 

segment (or pulse) lengths of turbulent and 
potential signal, respectively, in inter- 
mittent signal 
time 

velociU" along jr, ?/, and 0 , respective!}^ 

mean velocity on axis of jet or wake 
mean velocity in free stream of boundary 
layer or wake 

skin friction velocity, V r^/p 
velocity fluctuation along y, and 0 , 
res])ectively 

average velocity of propagation of turbu- 
lence front yelative to fluid (perpendic- 
ular to its own plane) 

Cartesian coordinates (x is measured from 
beginning of working section in bound- 
ary-layer case) 

Cartesian coordinates alined locally with 
turbulence front 
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'^0 

Y 


F: -V- V 


ai)])aroiit origin of wake, boundary layer, 
or jet 

iiistaiilaneous y -location of turbulence 
fi’ont for boundary layer and wake 

instantaneous — location of turbulence 
front of plane wake 


V = 

{x-x,r^ 

Vll2 


a,/3 


7 
A 

8 


8 * 

e 





p 

^(Y-yy 


ill wake 

value of y at which mean velocity defect 
is half maxinunn 

radian frequency of wall oscillation in 
model of laminer siqierlayer 
intermittency factor, relative time spent 
by a fixed probe in turbulent fluid 
instantaneous vector velocity jump across 
laminar superlayer 

boundary-hiyer thickness, the value of y 
at which 

bo u nd ar v-1 ay er ( 1 isj)l a(‘cm en t t h ick n ess 
thickness of laminar superlaver 
model superlayer thickness for mean and 
tluctuating vorticity, resjiectiv'ely 
momentum thickness of boundary hwer 
Lagran^^ian length scale in flow direetiom 

ru 

transversal Eulerian microscah^ of turbu- 
lence 

Eagrangian time microscale 

microscale of Y(t) times ( 

Lagrangian length microscale, v'\t 
viscosity coefhcient 
kinematic viscosity coefficient, ju/p 
total vorticity components in x— and 2 — 
(1 i rec t ions , resp ec t i vely 
total vorticity, H = 

vorticity fluctuation components along x, ?/, 
and 2, resj)ectiv7'ly 
vorticity fluctuation, E = 
density 

standard deviation of {V — )')=Vi 


O' it — 

T 

To 

X 

CO 


stress tensor 

time interval 
skin friction stress 

rate of dissipation of turbulent energy per 
unit mass of fluid 

Kolmogoroff (minimum) length, x = 
autocorrelation function of )"i(0 
autocorrelation functions of a, ar, mid b 
autocorrelation function of trigger output 
total vorticity vector 
vorticity fluctuation vector 


()pei*ators; 

( ) average 

( ) ‘‘short’’ time or s])ace average 


( )' = a ( y root mean square 

( hypothetical variable equal to actual varia- 

ble in turbulent fluid only and obtained 
by deleting potential fluid ]iart of an 
intermittent oscillogram 

EXPERIMENTAL EQUIPMENT AND PROCEDURES 

AERODYNAMIC EQUIPMENT 

ddie wind tunnel (fig. 2 ) is an open-return NEL type with 
a 2- by 2-foot working section and a free-stream tuidnileiice 
^evel at entrance of u'/7 '=0.05 ])ercent and ^;'/L =0.06 per- 
cent at a iin^an velocity of E =26 feet ])er second. 



Ficuue 2. — Schematic diagram of wind tunnel. 


In order to have a reasonably thick turbulent boundary 
layer in the relatively short working section, a wall was used 
as a working surface, and it was roughened by corrugated 
jiaper starting from the beginning of the contraction. The 
corrugations, set jierpendicular to the flow, were roughly 
sinusoidal, with about K-incli wave length and ^^4-iiich ampli- 
tude (half height). 

ddie exteJit of two-dimensionality in the boundary-layer 
flow was checked by mean velocity profiles at several stations 
across the 2~foot width of the working surface, at the farthest 
downstream station, j = 102 inches. The uniform zone was 
18 inches wide, with a boundarv-layer thickness of 
inches from wall to free-stream velocity and, estimating from 
reference 4, the transversal Eulerian scale was about 0.5 inch. 

The boundary-layei- measurements were all made at a 
free-stream velocity of 37 feet per second. The static pres- 
sure was very nearly constant along the working section 
(fig. 3). From comparisons with earlier work on this type 
of flow (ref. 11), it appears that the flow state is such as to 
have a fully rough wall condition. 

The round-jet unit is sketched schematically in figure 4. 
The orifice diameter was }i inch and it was run at an exit 
velocity of 300 feet per second. 



Firu'RE 3. — Static-])ressiir(‘ distribution along wind-tnnnel working 
section, (jo, dynamic pressure in free stream at x=d). 
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HOT-WIRE SET 

Most of the measurements reported here were made with 
the hot-wire anemometer as sensing element. The basic 
amplifier and compensation unit, constructed by Mr. C. L. 
Thiele, is described in reference 12. T\\o oscillograms were 
taken with a General Radio Type 761 camera photographing 
blue cathode-ray tubes. Measurements of the statistical 
distribution of lengths of turbulent bursts were made by 
scaling directly from the recorded oscillograms. 

The power spectra were measured with a Hewlett-Packard 
Type 300A wave analyzer, followed by a vacuum thermo- 
couple. The strongly fluctuating output was averaged by 
integrating with a fluxmeter and bucking circuit as illustrated 
in reference 12. 

The hot-wires used were either 0.00010 inch platinum or 
0.00015 inch tungsten, with lengths of about 1.5 millimeters 
for the i^-meters and 2 millimeters for the X-meters used to 
measure v', w\ and uv. No correction was applied for finite 
wire length. 


MEASUREMENT OF INTERMITTENCY 

Following Townsend (ref. 6) the intermittency 7 is defined 
as the fractional time spent by the (fixed) probe in turbulent 
fluid. Townsend has measured y in two ways: (a) from the 
‘^flattening factor’^ (or “kurtosis’O of the probability density 
of the intermittent signal (ref. 6); (b) from the mean-square 
output of an on-off signal triggered by passing the intermit- 
tent signal through a gate (ref. 10). The method used here 
is a development of (b), the relative “on time^^ being meas- 
ured by counting a high-f requeue}^ pulse signal as modulated 
by the on-off signal. This should give more accurate results 
at low values of 7. 

The overall block diagram is given in figure 5. Figure 6 is 
a further breakdown of the manipulative details, with a 
schematic diagram of a hypothetical signal as modified by 
passage through the various blocks. The actual circuit 
of this is given in figure 7. It is clear from figure 6 that the 
number of pulses counted for a given input signal will be 
a monotonically increasing function of discriminator setting. 
One would like to find a wide range of discriminator settings 
over which the count rate, for a given input signal, would 
be unchanged. Unfortunately, there is no such indic^ition 

Compensated 
input from hot-wire 
circuit 



Amp 

lifier 


Differentioting circuit 


Amp 

lifier 



Intermittency circuit 


Counter 

Figure 5. — Overall block diagram of intermittency-measuring 
arrangement. 
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Figure (>. — I4et ailed block diagram of intermit tency-nieasiiring device. 


computod from turbulonce level and microscale measure- 
ments. Estimates of the parasitic sensitivities, especially 
to the three (*omj)onents of turhulent velocity, were made h}^ 
measuring; the stead3^-state }"aw and speed sensitivities in a 
low-turhulence stream. These were found to be negligibly 
small for the particular meter used in getting the data. 

Xo correction has been made for finite wire length (the 
lengths were about 1 millimeter), and no correction has been 
made for the nonzero ratio of wire spacing to turbulence 
microscale, a characteristic giving parasitic sensitivity^ to the 
second derivatives of velocity fluctuations. 

MEAN VELOCITY PROFILES 

A flattened no. 20 hvpodermic needle was used as total- 
head lube in the measurement of the mean velocity profiles * 

from wliich boundarv-laver and jet thicknesses were de- | 

termined. 

Although exact wall location is ])robabl\^ a meaningless ^ 
concept for i*ough-wall boundarv-laver flows, the choice of | 
such a i-eferencc value of y is convenient for presentation of 
data in familiar coordinates. Therefore, a y =0 reference was 1 
chosen bv extrapolation to zero of the mean velocity^ profiles | 
from a region outside tlie boundarv tangent to the corruga- 
tion ])eaks. In order to minimize scatter near the ‘Svall,^^ ! 
all total-head traverses were made at the same phase position 
in the corrugation peak. A slight cutout on the downstream 
side in each case jicrmitted the total-head tube to go com- ^ 
])letely' into the boundarv. 

Since the exact details of mean velocity^ profile shape were i 
not of ])rirnarv concern in this investigation, no coi’rection | 

for the effect of turbulence has been applied to the total-head j 

tube data. | 

MEASUREMENTS 


of a ^k*orrect’^ setting for the discriminator, ])ossiblv liecause 
of the lag introduced in the necessary smoothing ])rocess. 
A typical illustration is given in figure 8. 

In practice, the discriminator level was set for each signal 
bv visual observation on a dual-lieam oscilloscope of simul- 
taneous traces of the difrerentiated hot-wire signal and the 
corresponding trigger output (e. g., fig. 9). The settings of 
the noise clipper and of the smoothing-filter time constant 
were chosen bv visual comparison at the beginning of the 
sequence of tests and kept fixed for the entire investigation. 

The intermit tencv circuit was desgined and built by Mr. 
Donald S. Johnson. 

VORTICITY FLUCTUATIONS 

fldie pyramidal configuration of four hot-wires connected 
in a Wheatstone bi*idge responding pi*imarilv to the vorticity^ 
fluctuation component along the flow direction is due to 
Kovasznay (ref. 18). Figure 10 is an isometric sketch and 
a wiring diagram. Some of the jiertinent details are given 
in reference 14. 

Calibration of sensitivity^ to vorticit.v has been tried by^ 
spinning the meter about its axis (ref. 13) in a uniform flow, 
but for the measurements jiresented here an indirect method 
was used: The readings in a decaying isotropic turbulence 
were compared with the values of vorticity' fluctuation level 


MEAN VELOCITY FIELDS 

Rough- wall boundary layer. Mean velocity' ])rotiles as 
determined from total-head tube measurements are plotted 
in dimensionless form in figure 1 1. Ofliere is reasonably' close : 
similarity'. Of course, exact similarity is not to be expected ' 
since boundary'-lay'cr Revnolds number varies considerably' , 
with X and effective roughness varies slightly'. 

Tlie moment um-tliickness distribution 

is given in figure 12. The solid line is a simple power law i 
drawn from the ajiiiarent origin x =Xo. The similarity shown 
in figure 1 1 is close enough so that the displacement tliickness ‘ 

and the total thickness 6 (j'), the value of y at which / '= f/„, | 

are assumed jiroportional to 6{x) for the purpose of later 
figures. The 6{x) values are assumed to be more reliable ' 
than 5* because equation (2) deenqihasizes the relatively 
uncertain region near the wall. The values of 5 *(j*) and b{x) \ 
are then given by 6{x) times the average values of o*/^ and ^ 
bjS. These values are 1.47 and 7.2, respectively'. 
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Fir.URE 7. — Circuit of interinittency-measuriiig device. Capacitances are in microfarads unless otherwise noted; resistances, in ohms. T*, 

test count rate. 


r 


2 


O. I 
Q. ' 
< 


J I L_ 


Optimum discriminator setting 
° [(visual inspection oscillograph) 

t O 


I I I 1 1_ 


Discriminator setting 


Fir.URE 8. — Intermittency reading versus discriminator setting for a 
typical probe position. 



(a) Discriminator set too high (signal is du/dt). 

(b) Good discriminator setting (signal is r(0)- 

Figure 9. — Oscillograms of hot-wire signal and trigger outj)ut. 
357386—56 2 




Figure 10. — Vorticity meter. 


The skin friction coefficient 


Cf = 



( 4 ) 


computed from d(x) by the Von Karmdn integral relation 


is included in figure 12. 



(fi) 


Round jet. — ^Figures 13 and 14 present data for the round 
jet corresponding to the data for the boundaiy layer. The 
tail depression is, of course, due to the directional sensitivity 
of the total-head tube; at the jet edge, the mean velocity 
is chiefiy radially inward. These measurements agree with 
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Figure 11. — Mean velocity profiles at several stations in rough-wall 

boundary layer. 



Figure 12. — Momentum thickness and skin friction coefficient in 
rough-wall boundary layer. 

the results of referenees 2 and 12 on velocity jirofile and 
linearity of jet inomentiim sj)read with x. However, the 
angle of spread is slightly greater than that in reference 2, 
being 10.8° total angle for the iialf-velocity cone as against 
9.5° in tlie earlier work 

Possible factors in this difference are tlie following: 

(a) Different orifice boundary conditions: In reference 2 
the jet emerged from a plane wall about 25 orifice diameters 
in width; here there was no wall. 




Figure 14. — Half-velocity radius of jet. 


(b) Different Reynolds numbers: At the same xjro this jet 
field has a Reynolds number five times bigger than the one 
in reference 2. 

(c) Different measuring insti*ument: The mean velocity 
])rofilesof refenuice 2 were measured with a hot-wire anemom- 
eter while these were measured with a total-head tube. 
In neither cas(‘ were the data corrected for the (diflVrent) 
effect of turbulence on ajiparent mean velocity. 

However, the difference is of no interest here since the 
])rincii)al concei’ii is a comparison of the relative behaviors 
of overall mean flow field and irregular turbulence front. 

CHARACTER OF FLUCTUATIONS 

The intermittent character of the outer part of the turbu- 
lent boundary layer is indicated by ty])ical oscillograms. 
Figure 9 includes u{f) and duldt, while figure 15 includes 
u(t) and $(t). Obviously there are still appreciable velocity 
fluctuations in the nonturbulent parts of the flow, ddiese 
are of relatively low frequency. The typical time record of 
vorticity fluctuations ^(f) indicates that the nonturlnilent 
parts are irrotational, since the order of magnitude (d the 
low-frequency fluctuations visible between turbulent seg- 
ments can be accounted for by parasitic sensitivity in this 
])articular vorticity meter. 





Figure 15. — Oscillograms of vorticity fluctuations and of longitudinal 
velocity fluctuations in intermittent zone. 

A definite property of the u{f) oscillograms is one-sidedness 
of the turbulent bursts. This result shows that on the 
average the bulges of turbulent fluid are moving more 
slowly than the nonturbulent fluid passing by the same 
lateral y-position in the boundary layei. This is not sur- 
prising, since such turbulent bulges must largely originate 
from further in toward the fully turbulent regioip which is a 
region of lower mean velocity in the l)oundarv-layer case. 

This one-sid('dness is shai’iiened up a bit b}" the fact that 
(as will ])e proved later) the irrotationally fluctuating fluid 
must be traveling at the same mean velocity as tlu‘ free 
stream. 
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The qualitative description of tlie turhuleiiee propagation 
phenomenon given in the ^Tntroduction’’ requires that it 
aetually takes place through a (presumably thin) viscous 
shear layer plastered all over the boundary. In fact, this 
^daminar superlayer’’ is the boundary between turbulent 
and nonturbulent fluid. Inspection of the oscillograms 
reveals no clear similarity among all the beginnings and ends 
of the turbulent bursts, but this is not a contradiction of the 
pln^sical picture. Any such tendency must be completely 
masked by the randomness of velocity gradients (and hence 
the shears) in the laminar superlayer. Furthermore, the 
boundary itself is an irregularly wrinkled surface in three 
dimensions so that the relative orientations of hot-wire 
and boundary at the moments of immersion and witlidrawal 
are also random. 

TURBULENCE LEVELS 

Turbulence-level distributions for the three velocity 
components u. jU ^ lil at the boundary-layer 

station studied in detail (x=\02 inches) are plotted in figure 
16, with the corresponding mean velocity profile included 
for reference. Clearly the velocity fluctuations due to the 
presence of the boimdaiy layer extend far outside the region 
conventionally identified as the boundary layer. 



X, i n. 

Figure 17. — Boundary-layer turbulence level as a function of x for 
corresponding y-positions 

the squared fluctuation intensities and {w'Y. 

The shear correlation coefficient Ruv =^wju'v' becomes quite 
uncertain in the outer part of the boundary layer because 
the measurement then involves the taking of small differ- 
ences between relatively large uncertain readings. 


0 

Figure 16. — Turbulonce-levol distributions at a;=l()2 inches in 
boundary layer. 

Variation of turbulence level in the r-direction, for cor- 
responding locations in the boundary layer, is indicated by 
v'jU versus x at several fixed values of yjb (fig. 17). Since 

the Prandtl friction velocity ^ 7 = -^/“ is probably the basic 

reference quantity with the dimensions of J^ength/Time in 
a solid-wall shear flow, one expects that, for corresponding 
positions in the boundary layer, u' , r', and w' ccl\. In 
turbulent pipe flow, Laufer (ref. 15) finds that v'jVj versus 
radius is independent of Ke\Tiolds number except in the 
vicinity of the wall. This suggests that, in the boundary 
layer, v'jUr versus yjb may be constant away from the wall. 
Figure 17 shows at least no clear-cut contradiction with 
this hypothesis, within the ovei’all expei-imental uncci*tainty. 

TURBULENT SHEAR STRESS 

The turbulent shear stress distribution —puv{y), at 
x = 102 inches, is presented in dimensionless form in figure 
18 and shows the same behavior as in the smooth-wall cases 
(refs. 4 and 16), approaching zero apjireciably faster than 



VORTICITY FLUCTUATION LEVEL 

ddie measured distribution of root-mean-square vorticity 
fluctuation (^', the r-component onh^) across the boundary 
layer at a: = 102 inches is given in figure 18. The instrument 
was by chance sufficiently symmetrical that, within the 
purposes of this investigation, no correction for parasitic 
sensitivities was necessary. 


.003 

002 

(R. 

' 0 ^ 

001 
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Figure 18. — Vorticity fluctuation and turbulent shear-stress distribu- 
tions in boundary layer at x=102 inches. 

INTERMITTENCY 

Boundary layer. — The transversal distributions of inter- 
mittency 7(y) at several x-stations in the boundaiy layer (typi- 
cal comparison with V\U ^ in fig. 19) show good similarity 
wlien y is normalized with o-(r), the square root of the 
second moment of dy/d?/ with }^-origin chosen so that by/dy 
has zero first moment (fig. 20). As pointed out in the ^Tn- 
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ti'oductioii,” dy/dy is tlio prohahility dcnsitA’- of Y{t), tlic 
instantaneous y-position of tlie front between turbulent and 
nonturbident (bad, at a fixed x. 'Pben 

crW=[(F^]‘"'^[ j'jy-Ff g d{y-Y)Y m 

whore ^7/5// is written as a funetion of {y — Y). Therefore, 
(j{x) is a suitable measure of the width of tlie intermit tent 
zone, that is, of the wrinkle am])litiide of the turbulence 
front. 




cr 


Fioire 20. — Iiitermittency distributions for several x-stations in 
boundary layer. 

Another important statistical measure of the turbulence 
front is its average location, 

vSince d 7 /d// turned out to be symmetrical and, in fact, 
virtuall}" Gaussian within the experimental precision (see 
section “Probability Density of Y{ty^), the determination 
of a and Y was considerably simplified. Both and 

F(«^) are given in figure 21 and 8(x) is included for compari- 
son. The logarithmic ])lot was used to estimate exj)onents 
in ])ower-Iaw approximations for the three quantities. 

"Idle ])ower-law fitting has been done with the best common 
origin for the three sets of points in order to simplify the 
comparison concept. 

Round jet. -Intermittency data for the round jet cor- 
responding to the data for the boundary layer are given in 
figures 22, 23, and 24. 




(a) Function of x—Xq. 

(h) Function of x. 

Fkujre 21. — .\in{)litud(‘ and average position of turbulence front in 
boundary layer as functions of x and x — Xo. 




r-R 

a~ 

F luuKU 23. — Intorinittency distributions for several x-stations in round 

jet. 
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(a) Amplitude of turbulence front. 

(b) Average position of turbulence front. 

Figure 24. — .\mplitude and average position of turbulence front in 
round jet as functions of x/2ro. 

Townsend’s plane wake. — For conveuieiit coinparison 
Townseiurs last published data (ref. 10) for the plane wake 
have been put into a form eorresponding; to that of the other 
data (figs. 25, 26, and 27). 

Since, however, only the points for x\d =800 and 950 are 
in the fully developed wake, no attempt has been made to 
determine separate power laws from his data. Instead, 



Figure 25. — Typical intermittency distribution across a plane wake. 
(Data from ref. 10.) 



Figure 20. — Intormittoncy distributions for several a;-stations in 
plane wake. (Data from ref. 10.) 



Figure 27. — Amplitude and average position of turbulence front in 
j)lane wake as functions of x. (Data from ref. 10.) 


])arabolas have been drawn with his choice of apparent 
origin siinpl}" to show that his results are not in contradiction 
with the parabolic and }^(j) (predicted theoretically 

in a later section). 

STATISTICAL ANALYSIS OF ON-OFF INTERMITTENCY SIGNAL 
(Output of Schmitt Trigger) 

As sketched in figure 6, one stage in the electrical signal 
manipulation sequence is a two-valued (on-off) random 
function. These flat-toj) pulses have duration equal to the 
time spent b}^ the hot-wire in turbulent fluid and spacing 
equal to the time spent in nonturbulent fluid. 

Two basic statistical characteristics of sucli a random 
on-off signal are (a) its power spectrum and (b) the prob- 
ability densities of its top lengths and its bottom lengths. 
E.xcept in special cases, no one has yet deduced a relation 
between these two functions (see section ^^Probabilit}^ 
Density of Pulse Lengths'')- 

Since the jumps in this signal are generated by the random 
occurrence of a particular amplitude of a more general sta- 
tionary random variable, that is, F(0, its properties give 
some information on the ])roperties of Y{t). For example, 
the probability densities of top and bottom lengths indicate 
the statistical distribution of wave lengths of the turbulence 
front, though less directly than the way in which dy/dy gives 
the statistical distribution of amplitudes. A detailed dis- 
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(Tissioii follows in the section ‘‘Statistical Description of 
Turbulence Front 

The ])ower spectriun of the on-otf signal must he related 
to that of the total hot-wire signal, though not in any simple 
fashion. As will he pointed out later, considering the total 
signal as (‘ontinuous turhulence modulated hy tliis on-otf 
signal, it appears that carrier and modulation must he 
statistically indejxMulent for the power spectra to comhine 
simply. 

Figure 28 is a series of j)ower spectra F///) of the Schmitt 
trigger output at various values of 7 for t = 102 inches. 
Statistical symmetiy of F(/) (indicated hy the api)i‘oximate 
symmetry of ^7/5?/) i*e(iuii*es that Fc{n) for* intPi-mittency 
7 = 7i he ecpial to /A(/0 for inteiMnittency 7=^1 —71. Figure 
29 gives the pi'ohahility densities of tops and of bottoms at 
the same hot-wii’e locations, ddiese wei*e obtained by dii*er*t 
measui-ement of oscillogi’ajihic I’ecoi’ds. 

The solid line in figui*e 28 is the theoi*etical power siiectrum 
for a random flat-top signal whose jumps have a Poisson 
distribution in time (see section “Powcu* Spectrum of Sdimitt 
dhlggei’ Output,’’ especially ec[. (88)). 



Figure 28. — Power spectra of Schmitt trigj^er outputs for three differ- 
ent intermittencies in boundary layer at x=102 inches. 

Fc{n)_ 1 

F,(0) l+2.49X10'-‘/P 

THEORETICAL EXISTENCE OF TURBULENCE FRONT 
Although the relatively shar]) front between turbulent 
and nonturbulent fluid has been well established exjieri- 
mentally, this a])])arently ubiquitous phenomenon must still 
be explained and exploi’ed analytically. The oscillograiihic 
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Figure 29. — Probability densities of segment lengths of intermittent 
signal for three different intermittencies in boundary layer at x=\02 
inches. 


recoixls indicate that it is likely to be a boundary between 
rotational and iri-otational motion. The theor’etical dis- 
cussion will th(*refor-e aim first at heuristic dernonstr’ation of 
this concept by showing in this context the known fact- that 
turbulent stretching of the vortex lines in a local vorticity 
gr-adient tends to steepen the gradient (leading, of course, in 
the limit to zero vorticity on one side). 

Succeeding sections will discuss some of the ramifications 
of this physical ])icture, in pr-eparation for the more detailed 
analyses which follow. The degree of agreement between the 
])redictions of these analyses and actual experimental rvsults 
will provide further indication of the validity of the hypoth- 
esis that the nontui-bulent field is aidually iiTotational. 


STEEPENINer OF A VORTICITY GRADIENT WITH LOCAL PRODUCTION OF 

VORTICITY 

Since the distinction to be made here between tuiFulent 
and nonturbulent flow is on the basis of ])re'sence or absence, 
respectively, of random vorticity fluctuations, the boundary 
])henomena must obviously be studied in terms of vorticity 
as a ])i’incii)al (haract eristic variable. 

The vector form of the vorticity equation for Ihi-ee- 
dirnensional inuornpiTssible viscous How is 




( 8 ) 
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wliero ^ = Stokes derivative (following a 

fluid element), Q is total vortieity vector, 0 is total velocity 
vector and v is kinematic viscosity. 

Introducing a Reynolds type restriction: 

Rfe j 0 = R (£) + ^(Zy 0 w = 0 

Q(r,t) = Q{r) + q(r,t) £= 0 

equation (8) yields the mean-vorticity equation: 


vorticity-dominated phenomena must he associated with the 
fine structure of the turbulence, especially for high values of 
turbulence ReATiolds number = where u' is root- 

mean-squa]-e velocity fluctuation in the ^-direction and X is 
the Eulerian microscale. For large enough values of R\ 
there should exist a time long compared with that charac- 
terizing the main body of vortieity fluctuations but short 
compared with tliat characterizing the largest scale velocity 
fluctuations, which dominate the convective properties of 
the turbulence. For example, one can expect 


(^.V)i2-h((/.V)w=(^-V)y+(w*V)r/+rV“^ (9) 

Subtracting equation (9) from equation (8) leaves the 
equation for vortieity fluctuation: 

(fi-V)(/+(w-V)^+(w-V)r/-(co-V)5 (10) 

The scalar product of w with equation (10) gives the 
equation for instantaneous vortieity intensity: 

I K£-v)2] [(^1 =co- [( 2 - v )£] + 

w-[(-wV)f^^] + w-[(cj*A)(/] — wd(^ii‘V)^] + rco-(V“w) (1 1) 

In Cai-tesian tensor notation, but keeping vortieity as a 
vector instead of an antisymmetric second-rank tensor. 


2 Dt 




5H, 

dxf 





dr, 

dxj~^ 




dXj 




(11 a) 


where a repeated index indicates summation and co‘^=^k^k- 
The averaged equation is 


1 77 . 1 d . I P — “ y 1 

9 ^ ) + ^i% —"At :^—+ 


dXj 2 dxj 


dXj 

— dr, 


dxj 


dui 




( 12 ) 


It was Ta\dor (ref. 17) who first identified the 

rate of production of vortieity fluctuations by the random 
stretching of vortex lines. It is largely the absence of this 
effect that makes fully two dimensional motion trivial in the 
problem of fully developed turbulence. 

To demonstrate the tendency of a vortieity gradient to 
steepen in the presence of this vortieity production effect, 
consider the simplified form of equation 1 1 la) for a flow with 
no mean velocity or vorticit}': 


do)2 .1 ^0)2 / d^,\_ 

bt "^2 dxj dxj) 


1 
9 




(13) 


It appears that no conclusion can be reached without further 
restriction. Since the vortieity spectrum varies like k^F{k)^ 
the running second moment of the velocity spectrum, 


where is the root-mean-square a:-component of vortieity 
fluctuation and X, is the Lagrangian time microscale (ref. 1). 
Introducing \ = v'\t (ref. 18) and, with local isotropyq the 

isoti’opic relation equation (14) becomes 

X 


0.45 



(15) 


or, in terms of A\, 
(ref. 18) 


the large R\ approximation for \j\ gives 



(16) 


For flows with equation (16) valid, equation (13) could be 
averaged over a time long enough to average vortieity 
phenomena but short for convective velocity phenomena: 


1 1 ^^2 

2 U "^2 dX) 




(17) 


where ( ) « ( ) for the fine-structure variables. 

The velocity derivative has characteristic time like that 
of vortieity. 

Since the objective is to show the steepening of the 
gradient in the absence of v, omit the last term and write 


Therefore, 


1 d^ 

dXi 


1 D (b7^\ 

2Dt\dyJ^dyV'^^ dxj 


(18) 


(19) 


Taylor (ref. 17) has shown that ^^>0 for isotropic 

Ox j 

turbulence. In fact, since this inequality just expresses the 
general tendency for fluid lines to lengthen in a turbulent 

flow, it seems clear that ^^>0 in any turbulence. Then, 
if monotonic with c/, it follows from equation (19) 

OXj 



has the same sign as do:-/dy, which means a 


steepening of this gradient. 
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Tlu‘ ])hysical roasoii lor tlie steo])ening tondoiicy is, of 
course, just the fact that tlie I’ate of production of new 
vorticity l)y line stretching is ])roportionaI to the vorticity 
already ])resent at any point in tlie fluid. Hence the higher 
vorticity regions ex])erience a greater i-ate of increase of 
vorticity than the lower vorticity regions, that is, the 
gradients tend to stee])en up, limited finally by viscous 
ditVusion and dissijiation. Of coursi‘, the gradient of con- 
cern here is that in the zone between fully turbulent fluid and 
nonturbulent fluid. A steepening of this gradient means a 
tendency toward a relatively shar]) surface of demairation 
between tlu' two states. The alcove discussion does not 
treat the (piestion ol the ecpiilibrium thickness eof the laminar 
superlayer that results; this will be estimated later. Of 
course, for the turbulence front to be sharp as observed 
experimentally, it must be shown that e^<5. 

Although the analysis is valid only for ewtremely high 
values of /i\, far higher, in fact, than those that occur in the 
expiu'iments reported here, there appears to be no reason for 
the situation to change c[ualitatively at knver values of /A, 
as long as nonlinear effects in the Navier-Stokes ecfuations 
remain ini])ortant, for example, 

LAMINAR SUPERLAYER 

Vorticity can l)e transmit t('d to an irrotational flow^ only 
through th(^ tangential forces due to viscosity; it cannot be 
transmitted to the irrotational flow^ by macroscopic Reyn- 
olds type shear forces. It therefore follows that the 
instantaneous border zone lying between turbulent fluid and 
irrotational fluid must be a region in which viscous forces 
play a central role, in s])ite of the presence of velocity 
fluctuations which dominate the gross momentum transfer 
of the turbulent field. This border zone may be termed the 
laminar su]K‘rlayer and is exactly what is also referred to in 
this report as the turbulence front, although the latter 
designation inpdies emjdiasis on its overall beha\ior rather 
than its detailed structure. 

This laminar superlayer differs in function from the wull- 
known laminar sublayer at the smooth solid boundary of a 
channel, l)ipo, or boundary -layer flow\ The sublayer is a 
relativ('ly fixed region in wdiich mean flow* momentum is 
trans])orted primarily by a net mean viscous (laminar) 
shear force. It transmits \jii\e mean vorticity (being a 
zone of roughly constant il(y)) and it remains ^hittacluar' 
more or less to the same fluid particles. On the other 
hand, the superlaAVu* is a (convectiv(‘ly) randomly mov- 
ing layer of fluid which probably transpoids relatively 
small amounts of mean momentum and vorticity by viscous 
shear forces; its distinguishing function is transport of 
vorticity fluctuations and mean vorticity, wdien present, 
into what w^as previously an irrotational field, and in so 
doing it continuously propagates (relative^ to local fluid) 
normal to its local ^^plane.’^ 

Sketches (a) and (b) illustrate the concept of the super- 
layer as a very narrow’ zone in wdiich the vorticity fluctuation 
level and the total shear (if any) drop from values charac- 
teristic of fully turbulent flow’ to practicalh^ zero. 


Fluctuating potential flow 


^^Laminar supe'layer 



U{y) 


Turbulent flow 
Sketch (a). 



While the instantaneous local viscous shear force in"[a 
laminar sublayer is pnalominantly in the direction of the 
mean shear force, that in the superhwer must have a much 
higher fluctuation level, often reversing its direction, for 
exam])le. In fact, in a flow’ field w’ith constant mean velocity 
everywhere, the su])erlayer viscous shear force would have 
no mean valiu* at all. 

The discussion headed ‘‘Steepening of a Vorticity Gradient 
AVith Local Production of Ak’^rticity’^ is a justification 
(not a proof) of the experimental fact tluit the continuous 
fluid-line stretching due to the velocity fluctuations tends 
to stee])en up the laminar superlayer. This steejiening 
effect is reinforced by the projiagation and must, of course, 
be balanced out at some state by the diffusive action of 
viscosity, so that the sujic'rlayer must have some average 
thickness. Tiom the oscillograms, it appears that this 
c[uantity, e say, is very small. 

Some heuristic comments can be made about tiiis tliic km^ss. 
First of all, since tlu' layer is primai*ily a vorticity-propagating 
device, its thickness should be h'ss than a length chai*acter- 
izing vorticity fluctuations on the turbulent side of the 
boundary, foi* example, the dissipation scale X (Taylor^s 
microscale'). In fact, as a characteristic viscous shear 
length, it might be expectc'd to be the same order as 

Kolmogoroff’s minimum length ( — 1 ? whei’e is the rate 

of dissipation of turbulent energy per unit mass of fluifl. 

A second intuitive s])ecification is that, as a violently 
disturbed free laminar shear layer, its (‘haracteristic Reynolds 
number should be on the order of the low’er critical Reynolds 
number for free laminar shear layers. A possible choice of 
characteristic instantaneous Reynolds nundier would be 
that based on thickness and instantaneous tangential 
velocity difference A across the superla 3 ^er. AAlien the 
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instantaDoous viscous shear force, S^/jL = por unit area in 

(lie plane of (lie front, has a nonzero average (i. e., a pre- 
ferred direction), a plausible average Keynolds nuinber 
might be 

( 20 ) 

wliere A is the magnitude of A. Of course, A and e are 
doubtless negatively correlated, but inclusion of such a 
refinement would be inconsistent with the crude nature 
of the discussion. 

For turbulence fronts in which there is little or no mean 
velocity difference across the superlayer, the above definition 
is inapplicable and might be replaced by 

( 21 ) 

V 


A simple mathematical model of the laminar superlayer 
wdll be taken up as a separate section in the discussion of 
propagation velocit}" of (he turbulence front. 

The follow ing important inference can be made on the basis 
of the highly localized character of the laminar superlayer: 
Since no appreciable viscous effects extend be}’^ond this thin 
layer, and since only viscous effects can transmit vorticity, 
it follows that the mean velocity everywhere in the potential 
jmrt o f the flow must he constant and equal to that at ^^inflnityy 
This is a consequence of the fact that the mean vorticity is 
12=dL7d?/. 

This conclusion wdll be analytically emphasized in the 
followTiig section. It is in contrast wnth an assumption of 
Towmsend (ref. 10) that the nonturbulent fluid lying between 
bulges in the turbulence front ^fis constrained by pressure 
gradients to move at the same mean velocity^’ as the fluid 
in the adjacent turbulent bulges. 


again omitting the implications of Ae correlation. 

Since, how^ever, there still exists no analysis relating A 
to the properties of the turbulence, a third definition, re- 
placing i, 2 R« and including such properties, is preferable: 


In fact, this definition is not too different from the other two: 
A/e must be of the same order as the neighboring turbulent 
vorticity fluctuations. 

Given an order of magnitude of the Iowtt critical Reynolds 
number for free laminar shear flow' ])lus a measurement or 

estimate of oj' = iri the turbulence near the front, 

an estimate can be made for e. 

The only information available for estimating the desired 
Reynolds number is the partial analysis of Lessen (ref. 19), 
a small-perturbation anah'sis. Extrapolation of his neutral 
stability curve (a highly inaccurate process) suggests an 
estimate 

1</L<10 

The measured turbulent value of for a typical case 
(fig. 18) is about 400 per second, wdiich gives co'«700 per 
second, if there is approximately local isotropy. With 
r=0.15 square centimeter per second the estimate of super- 
layer thickness turns out to be 

0.015<e<0.05 centimeters 


IRROTATIONALITY AND REYNOLDS SHEAR 


In view' of the evidence that the fluctuations outside a 
turbulence front are irrotational, it is pertinent to take a 
look at the customary turbulent-flow equations (actually 
valid for anv stationary fluctuations) for the particular case 
of irrotational fluctuations. The hope is that some drastic 
simplification will appear. 

In Cartesian tensor notation, the Reynolds equation for 
steady mean motion is 


1 dP , d / X 

^ * dx, “ p ^ * bxj 


(23) 


Tlio last term is the turbulent apparent force vector or 
Revnolds vector. 


P"or irrotational fluctuations. 


duj_duj 
dxj dxi 

therefore, 

A (— .1 e>(^ _i dr/ 

2 dXi ~2bxi 


(24) 


which shows that the Revnolds force reduces to a normal 
force only, since it is expressible as the gradient of a scalar. 
The Keynolds equation can then be written 



This appears to be a reasonable order of magnitude since 
X«0.2 centimeter in this ])art of tlie flow . The Kolmogoroff 

lengtli roughly O.Oo centimeter. 

In concluding this section it should be mentioned that, 
although no systematic measurements of A have been .made, 
rough estimates from oscillograms in the intermittent zone 
of the boundary layer indicated the order of 0.05 to 0.10 
times t/„. This average velocity defect indicates the 
obvious fact that turbulent boundary bulges originate in a 
region of low er mean velocity and also represents the presence 
of vorticity and of locall}' laminar shear. 

357386—56 3 


and this form emphasizes the fact that, regardless of the mean 
velocity field, irrotational fluctuations yive no net apparent 
shear forces on a fluid element.'^ 

This does not necessarily mean, howTver, that the Reynolds 
shear force on a plane is zero or that the Reynolds stress 
tensor —pupTk has only leading diagonal terms. Also there 
may still be a continuous “production’’ of fluctuating kinetic 
energy, that is, a transfer from the mean motion kinetic 

energy (ulil-I^Vo). 


3 This fact was pointed out by Dr. F. H. Clauser, 
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Comparison of equations (23) and (25) for the ease of 
motion two-dimensional in tlie mean (three-dimensional 
irrotational fliietiiations) yields the relations for the Reynolds 
shear force comiionents 


duv Id,— ,— 

<);, -2 5 i 




duv 

dx 


1 d 

=2 57 /(^'+^- 


v^) 


(20) 


(27) 


Probably the highest intensity random irrotational tiue- 
t nations easily available in the laboratory are those in the 
“potential eone” of a round turbulent jet. These apparently 
get as high as u'lU^o percent (ref. 2). 

Ecpuition (28) also can be deduced for the special case of 
a constant mean velocity held with arbitrary tluctuations, 
provided only that the mean values are plane, that is, 


which may also be regarded as a pair of ditrerential equations 
relating the four nonzero components of the Reynolds stress 
tensor. 

dhvo provocative forms follow troni alternative combina- 
tions of equations (26) and (27): 


and 


d'nTv d‘^uv 
d/^ dif 


d^ 


dx dy 


(u“— r) 


dhiv ^d'-uv _ d~w' 
dx' dif dx dy 


(28) 


(29) 


Equation (28) gives the interesting conclusion that if —v^ 
is constant in either x or y the turbulent shear stress satisfies 
a homogeneous plane-wave equation with characteristics at 
±45° in the J?/-])lane. 

For the particular flows studied in this rejiort, the experi- 
mental results show that uu) ap])roa(*hes zero fastei* than u'v' 
as ?/ (or r) is increased. This seems to indicate that uv=^0 
in the potential field. However, insufficient coverage and 
accuracy of the data preclude the ])Ossibility of checking 
this through equation (28). Since ^ =u^ =0 for y=^, 

such a check would require that =v^ throughout the po- 
tential field. 

Parenthetically, viscous fluids with zero net shear force on 
a fluid element hut with nonzero shear stress are far from 
unknown mathematically: Any irrotational laminar flow of 
a viscous fluid is such a case (aside from the trivial case of 
f2 ^Constant). The requirement on a stress tensor 
that it produce only normal forces is that 


dai/c d G 

dxjc~dxi 


(30) 


where (7 is a scalar. 

The principal significance of equation (25) in the general 
problem under investigation is as follows: Assuming that 
the fluctuations on the free-stream side of the turbulence 
front are actually irrotational, as both measurements and 
heuristic reasoning indicate, the mean velocity there must 
be equal to that for y = «? . This verifies the physical infer- 
ence drawn in the previous section from the concept of the 
localized laminar superlayer. 

It appears paradoxical that the mean flow kinetic energy 
should be unchanged in a zone where there has appeared an 
appreciable kinetic energy in velocity fluctuations. How- 
ever, the latter can come from the turbulent part of the field 
through nonviscous effects, leaving mean flow kinetic energy 
in the potential zone unchanged. This would he consistent 
with the inference that i7F=0. 


THEORETICAL BEHAVIOR OF TURBULENCE FRONT 

As mentioned in the “Introduction,’^ two of the fluid 
mechanically pertinent characteristics of the relatively 
sharp boundary between turbulent and nonturbulent fluid 
are (a) its mean rate of increase of wrinkle amplitude in the 
downstream direction and (b) its mean velocity of propaga- 
tion transversely into the irrotational fluid. The following 
sections represent crude theoretical attempts to predict 
these two cliaracteristics in terms of the statistical properties 
of the fully turbulent fluid on one side of the boundary. 

WRlNKLINCi RATE 

In turbulent flows with greater than about 10, th<u*e is 
no reason to expect any ])articular chunk of fluid to return 
to the nonturbulent state once it has become turbulent. 
Therefore, the presence of turbulence in a small piece of 
fluid can be regar<led as an indelible tagging, somewhat like 
heat or a chemical contaminant. Were it not for the con- 
tinuous propagation of the turbulence front into new (luid, 
this front would always consist of the same fluid particles 
and would obviously be susceptible to a Lagrangian study 
in terms of Taylor's theory of diffusion by continuous move- 
ments (ref. 20), as has f)een applied to the wrinkling rate 
(identical to turbulent diffusion rate) of a very thin sheet 
of thermally tagged fluid in a turbulent flow (refs. 1 and 18). 

In fact, a uniform translational velocit}^ of the tagging 
attribute relative to the fluid does not render Taylor’s con- 
cepts invalid; it does, however, require a generalization of 
the analysis to a mixed Plulerian and Lagrangian treatment, 
though somewhat different from the relative dispersion case 
set up by Brier (ref. 21) and by Batchelor (ref. 22). 
Clearly in the limit of (e. g., wrinkling of a Mach 

wave propagating through low-s])eed turbulence) it recuces 
to a simj)le Eulerian diffusion problem, while in the limit of 
V*<^v' the pureU^ laigrangian analysis of Taylor applies. 

P^or the present problem it appears that neither of these 
limiting conditions holds, although the latter is closer. 
Consider the rough- wall boundary layer at j: = 102 inches 
as an example: 




v' dx 


where is taken in tlie fully turbulent zone adjacent to 
the intermittent zone. This formula is deduced in the 
section “Applications to Particular Turbulent Plows.” 

Since the basic jiroblem (diffusion of a front propagating 
through a homogeneous turbulence) has yet to be analyzed, 
the present phenomenon will be estimated as though 
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Subsequent approximations are eonsistently rough. 

Another peculiar property of the present problem is that 
the surface whose turbulent diffusion is of interest has 
turbulent flow on only one side; the thermally tagged sur- 
face used in conventional diffusion studies has the same 
kind of turbulence on both sides. However, the theoiy of 
diffusion by continuous movements is simply a kinematic 
anah'sis based on the presumably given velocity statistics 
of the fluid particles in the surface. If these are correcth" 
given, no further infoi*mation or restriction is necessary. 
Therefore, since the purpose of tliis section is to predict tlie 

form of a{x)=->i{Y-Yf in terms of the properties of the 
fully turbulent zone, the only additional assumption neces- 
sary is that the velocity fluctuations of the fluid particles in 
the front are i)roportional to those in the fully turbulent fluid 
near the front. 

ddie anaHsis of one-dimensional diffusion by continuous 
movements for a homogeneous field with no mean motion 
leads to 

^Jlr.{r)dr (31) 

where a is the standard deviation of the distance traveled due 
to turbulent convection and is the Lagrangian correlation 
(‘oefficient. For times long compared with that for vhich 
the familiar asymptotic form results: 

cr(0 = y'V^ (32) 

00 

where \ is the Lagrangian (time) scale. 

If a relativel}^ high uniform mean velocity in the ^-direction 
is introduced equations (dl ) and (32) can be inter- 

])reted approximately in s])atial terms since Ut^x for any 
particle (refs. 1, 18, and 23). Then 

(33) 

u 

where Al = ULi is an approximate longitudinal Lagrangian 
length scale. 

It has been pointed out in previous publications (refs. 18 
and 24) that the most concise representation in such a flow 
follows from introduction of a transversal Lagrangian length 
scale Ll^v'Li. Then 

o-(/) « -y /2 (^=^ LlX (34) 

whicli gives tlic dispersion (identical to surface wrinkle 
amplitude) at large distances from a fixed source of tagging, 
when v{t) following a fluid particle is a stationary random 
variable and U is constant. 

For the hypothetical case of the turbulence front bounding 
a turbulent motion homogeneous in the stream direction, this 
asymptotic form would pertain; the ‘‘source^^ lies indefinitely 
far upstream. However, in virtualH all turbulent flows of 
interest, the statistical properties of the motion vary with x. 
Consequently, application of equation (34) to these cases 
implies the further restriction that these o*-variations be slow, 
that is, that there be little change in an a:-interval comparable 
with A^. 


A particular example of the degree of validity of this 
restriction can be drawn from the case of decaying isotropic 
turbulence, where Lagrangian scales have actually been 
measured (ref. 18). At 43 mesh lengths behind a 1* inch- 
square mesh grid of }^-inch dowel, with U=25.8 feet per 
second, it is found that v'lU =2.0 percent, 17 inches, and 
(Iv'Idx corresponds to a change of about }i in v'lU over an 


X* interval equal to A^^. 


However, 


1 (ILl 
Ll dx 


Kl is only on the 


order of 0.03. 


In most shear flows, the v'lU changes will be slower than 
for this decaying isotropic turbulence while the Ll changes 
may be slightl^^ faster. In general, it can be anticij)ated that 
in the application of equation (34) to boundary layer, jet, 
and wake the requirement of slow ic-variations in turbulence 
properties will be satisfied at least as well as the previously 
mentioned restrictions for this Lagrangian treatment. These 
applications and comparison of computed values of (t{x) with 
experimental results will be presented further along, under 
the appropriate section headings. 


PROPAGATION VELOCITY BY DIMENSIONAL REASONING 

The average velocity of propagation of the laminar super- 
layer (or turbulence front) relative to the local fluid F* must 
be monotonic with the average magnitude of the instan- 
taneous (laminar) shear stress in the superlayer ^proportional 

to M However, the ratio is not a directly measurable 

quantity and must be replaced by something more tractable. 
As has been mentioned in the section ‘d^aminar Superlayer,’’ 
when there is no mean shear stress A/e must be of the same 
order as the vorticity fluctuations in the turbulent fluid near 
the front. Therefore F* should be monotonic in co'. Since 
this is a viscous phenomenon, it must also depend upon v. 

In fact, the inference that F*=F*(r, w') can be made on a 
much more direct and superficial level. Since the laminar 
superlayer is a device for the viscous propagation of vorticity 
fluctuations into an irrotational fluid (in the case of zero 
mean shear), the propagation velocity must depend at least 
on co' and on v. Furthermore, these alone are sufficient to 
pi*oduce a parameter with the dimensions of velocity. 

The only combination giving the appropriate dimensions 
gives, by inspection, 

F*oc ^ V 03 ' (35) 

for zero mean shear stress.^ Of course, F* is directed per- 
pendicular to the local tangent plane of the turbulence front. 
The effects of nonplanarity of the whole front will be noted 
later in this section. 

p]quation (35) would be expected to apply, for example, 
in the case of the boundary between a homogeneous tur- 
bulence and a nonturbulent fluid, with C constant over the 
entire flow field. 

At the free boundary of a turbulent shear flow it is to be 
expected that the shear force vector of the laminar super- 
layer will have a mean value which will also promote F*. 

« It should also be noted that the assumption that R«(where /?. = eC7>') has a “universal" 
average value corresponding to a lower critical Reynolds number coincides with the plausible 
dimensional hypothesis that eoc vWo'. 
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If the whole front were nearly flat, this mean value would 
he a function of y in the superlayer and also proportional to 
the mean shear in the turbulent fluid just inside the tur- 
bulence front, ^ varying from equality on the turbulent edge 
to zero on the free-stream side. Sketch (b) includes this 
concept in a coordinate system attaclied to the laminar 
superlayer. 

In this more general case, the physical picture suggests 
that V* depends upon the average magnitude of Uie total 
shear in the superlayer: where ^(0 = 8+.s*(/) is 

the shear force vector on a unit area on the turbulent edge 
of the superlayer. With Cartesian coordinate system Xu ?/i, 
and Zi fixed in and alined with the turbulence front (?/i per- 
pendicular to front), the fluctuation £ has only Xi- and 
^i-components. Then, with gross mean shear directed along 

Ji, 

Dimensional reasoning gives 



If the random slope of the turbulence front in the x-, ?/-, 
and ^-coordinates is small on the average, the X\-, ?/]-, and 
3i-system can be replaced by x, y, and z, and Sx^ is propor- 
tional to the mean shear stress in the turbulence. l^\irther- 

niore, with local isotropy in the turbulence, — As in 
the simpler case, these are proportional to Then 

equation (^16) can be written 

si+BAc^y 

where B is a numerical constant, probably of order unity. 
This reduces to equation (35) for a shear-free tijrbulence. 
For each particular type of turbulent shear flow, Sx can be 
taken proportional to some characteristic mean shear stress. 
Xo application of equation (36) or (37) is made later in 
this rejiort. 

Handling of the propagation problem in terms of a plane 
turbulence front implies that e is much smaller than the 
radii of curvature of the front. The degree of validity of 
this assumption is not easy to check directly from the sta- 
tistics of the turbulent fluid; it requires fairly detailed infor- 
mation on Y{x,f). However, the measurements on statisti- 
cal distribution of pulse lengths cmning out of the trigger 
circuit (fig. 29), transformed by from time to length, 
give indirect indication that the assumption is well satisfied. 

Conversely, since C* is normal j)ro])agation velocity of 
the front (especially in the case with zero mean shear), 
propagation with constant V* over the wliole front would 
tend to introduce a skewness into the jirobability density of 
y, as in sketch (c). This is the effect mentioned by Kaiiovitz 
(ref. 26) in accounting for the skew nature of the flame front 
as observed in a turbulent bunsen flame. 

Tlie highly symmetrical shape of byjdy (indicated by 
linearity in fig. 32) shows that this effect, if present, is 
negligible in the phenomenon considered here. 

* In fact, I\ II. Clausor proposes a propagation velocity, for the tiirb^ilent boundary layer, 
depending only on the mean shear stress in the turbulence: V*(xtoIpA (ref. 25). 




Potential 



Sketch (c). 


This negligibility is an indication that the radii of curva- 
ture of the front are large compared with the wave lengths. 
This means that V* is directed very nearly perpendicular to 
the }\x) surface. For two-dimensional flows in which the 
boundary-layer approximation applies, this in turn is nearl}' 
jiarallel to the j!^-])lane, that is, (IYldx<^\. Therefore, within 
a corresponding approximation, the surface area of the 
turbulence front on a two-dimensional flow is equal to its 
])rojection on the j:£-plane, and the average rate of conquest 
of ne\v fluid by the turbulent state is V*, in units of volume 
])er unit time ]>er unit area of contact. A similar concept 
holds for the axially symmetric flows. 

It is obvious that a turbulent shear flow can have similarity 
only if both (j{x) and Y{x) are proportional to the boundary- 
layer thickness 8{x) (wdiich, of course, must be proportional 
to any other characteristic thickness defined in terms of the 
mean velocity profile). Stated in different but related 
terms, the average rate of flow of turbulent fluid passing 
through any constant j-plane must be proportional to the 
rate of flow of boundary -layer fluid passing through the plane. 


MODEL OF LAMINAR SUPERLAYER 


Dimensional reasoning as employed in the preceding sec- 
tion, and in earlier ones, gives at best the functional forms 
of the laminar-superlayer characteristics in terms of the 
statistical properties of the turbulence with which it is 
associated. Fully quantitative results follow only through 
deductive analysis, that is, actual solution of an appropriate 
boundary-value problem. Since the actual problem appears 
to be too complex for full solution at present, a simple 
jiliysicomathematical model will be used with the expectation 
that the results, after interpretation in terms of pertinent 
variables in the actual problem, will give a proper order-of- 
magnitude relation among these variables. 

The model proposed is a generalization of the Stokes and 
Kayleigh jirobhmi of the infinite wall moving in its own 
plane (ref. 27). The first extension is the addition of a 
constant suction velocity r(<0) with, of course, wall 
])orosity. The velocity I" corresponds to propagation ve- 
locity of the turbulence front. The differential equations 
are thus 


dr 

bt^ dy bf 


(38a) 


dir dm 

di+ by~'' bi/ 


(38b) 


Since the U and W equations are independent they can be 
treated separately. In the absence of mean shear they are 
identical, and only one need be considered. 
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Since equations (38) are linear, the vortieit}" coinponents 
Z^dUldy and E = dHydy obey the same equations as the 
velocities: 


dt by " bi/ 


dU by bx/ 


(39a) 


(39b) 


'J'hesc eciiiations are to be solved witb boundary conditions 
Z(t,-)=H(<,co)=0 (40a) 


Z(<, 0)=Z<,+ C„ sin at 
2 (<, 0) = ^„ sin (8/ 


(40b) 

(40c 


Since E has onlj' a fluctuating part, and since linearit^y 
permits separation of tlie steady and fluctuating parts of Z, 
tlie problem becomes 


(l^Z 

-clxj-" dxf 


with Z(oo) = 0 and Z(0)=Z<„ and 


dy 


dif 


(41) 


(42) 


with f(/., oo) = 0 and f(/,0) = fo shi od. There is an identical 
boundary-value pi'oblem for 

To get closer equivalence to the fully three dimensional 
problem a purely mathematical extension can be made, cor- 
responding roughly to the pln^sical phenomenon of contin- 
uous vorticit}^- fluctuation production (by fluid-line stretch- 
ing) at a rate proportional to that already present. This 
is most simpl}^ done by adding a linear term to the ^ equa- 
tion, giving 


— + V — = V 

bt^ by bx/^ ^ 


(43) 


with a similar expression for Here the negative root has 
been chosen so that equation (45) reduces to Stokes^ solu- 
tion for H=/v = 0. 

From equations (44) and (45) it is desirable to extract an 
ox])ression for tlie thickness of tlie disturbed laj^er. A con- 
venient measure of thickness is simply the inverse of the 
coefricient of — ?/ in the exponentials of both solutions: 


V 


Cl 


V 


(46) 






1/2 


(47) 


Application of equations (46) and (47) to _tlie laminar- 
superla^^er problem requires identification of F, a, and K 
with measurable variables in the turbulent fluid near the 
superla}^: 

(1) the propagation velocit 3 L 

(2) the root-mean-square value of any one of the 

three orthogonal turbulent voilicitA' fluctuation compo- 
nents. In other words, root-mean-square vorticitA^ he 
regarded as a characteristic frequency of turbulence. For 
large values of A\, Iw local isotropy. 

with local isotrop}^, is a measure of the rate of 

fluid-line stretching. 

Substituted into equations (46) and (47), these give 





1 


(48) 




1/2 


(49) 


where K is like a constant average vortex-line stretching 
rate. 

Xo corresponding term is added to equation (41) because 
the plane form of the mean-vorticit}^ equation for turbulent 
flow (eq. (9)) shows no term identifiable as production of 
mean vorticity due to random turbulent stretching of 
vortex lines. 

The solution of equation (41) is 

Z=Z„ exp (y 7/) (T<0) (44) 

The solution of equation (43) is 


®-r)4-4-A{'t©W- 

Q+f}"y](F<0; cx,K>0) 


A simpler, more approximate form for €2 is attained after 
inspection of the experimental orders of magnitude of V* 
and For example, at the inner side of the intermittent 
zone, in the rough-wall boundaiy la^'er at x=\02 inches, 
V* ^ 1 .3 inches per second and —400 per second. Therefore, 
it turns out that equation (49) can be simplified hy liberal 
omplo^unent of chopped-off binominal expansions. The 
roughest (and simplest) resulting estimate is 

(50) 

Since the laminar superla^^er can be assumed to exist even 
in the absence of a mean vorticiU' field, it is reasonable to 
assume that the fluctuating part of the superlayer model is 
the more pertinent one. Then one ma^' take €2 to be e, giving 
the theoretical prediction 

(50a) 




( 45 ) 


as an order of magnitude. This is consistent with the earlier 
conjecture on the constancy and order of a possible Re^molds 
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number = It does not appear to he siiseeptihle to 

direct experimental verification, but, as mentioned earlier, is 
of the same order as the KolmogorofF (minimum) length 

(5,) 

which follows from ecpiation (50a) and the relation between 
X and X. 

Xo estimate of V* follows from ecpiation (50) and, insofar 
as a strictly fluctuating laminar superlayer is concerned, the 
dimensionally induced equation (35) remains as sole redic- 
tion of propagation velocity. 

However, equation (48) for the mean thickness gives 


With this representation, the total signal expressed in 
terms of the other quantities is 

a{t) = b{t}[aT{t) + A] — bar— Ay (54) 

since a = 0 by definition and b{aT+A) = baT-\-Ay. Whence, 
the mean-square value can be written as 

'^^=h'aA+2Ab''aT—{baTy^2AybaT+yA—y)A^ (55) 

The objective is to express as a function of a'^ and other 
necessary parameters. Obviously this is impossible without 
introducing some further restrictions, especially on the 
statistical relation between ariO and b{f). Therefore, 
assume 


(52) 

If there is a single layer, V^* should be the same for both fluc- 
tuating and average vorticity. If ei happened to be of the 
same order as €2, equations (51) and (52) would give V* ^ 
O(a^), but there seems to be insufhcient a ])riori basis to 
make this guess a formal part of tlie analysis. 

INFERENCE OF TURBULENCE PROPERTIES FROM INTERMITTENT SIGNAL 

Townsend (refs. 10 and 28) has suggested that it may be 
possible to compute the statistical properties of the turbu- 
lence inside the convex bulges of the turbulence front from a 
knowledge of the corresponding statistical properties of the 
full intermittent signal plus the intermit tency factor 7. His 
hypothesis is that, in effect. 



wlu're a{f) is a random property of the flow and Ot is the same 
property but confined to the turbulent par^s o£th^total sig- 


nal. 


For example. 


Townsend refers 


to as 

7 


^‘the 


mean turbulent intensity within the jets” (identical to 
hulges). 

Actually, the applicability of eciuation (53) is contingent 
u])on very definite restrictions. For discussion jiurposes, 
suppose tiiat a(t) is the complete signal a_nd b{t) is the inter- 
mittent (0 or 1) signal. Obviously, b^b~=y, the intermit- 
tency factor. Also aAt') is a hypothetical signal whose ])hys- 
ical nature is the same as a(t) but ajiplies to turbulent fluid 
only. If Gt is chosen to have a mean value of zero, then it is 
necessary to introduce a constant quantity .1 which is the 
distance between the zero line of ar(.t) and the signal level 
cori*esponding to b(f) at zero. 

Sketch (d) illustrates the definitions. Implicit in this for- 
mulation and sketch is the restriction that a (or a^) is a physi- 
cal property which is zero in the ])otential flow region. 



(a) 607^=0 (whence b^ar—O, since 6^=6 ) 

y ( 56 ) 

(b) b~aT~ = b^XaT‘^=yaT^ J 

A sufficient but not necessary condition for these two is 
that b(t) and ar(t) be statistically independent. 

With restrictions (5b), equation (55) reduces to 

a^=yar-^Ay(l—y)A'^ (57) 

and the turbulence property ^ can be computed from the 
corresponding total-signal property plus measurements of 
7 and A. 

For some physical variables a(f) it will turn out that 
A =0 and then equation (57) reduces to eciuation (53). 

For the quantity to have any simple interpretation it 
must of course be assumed that the physical variable it 
represents is a homogeneous random variable in the turbulent 
fluid. 

Summarizing the conditions necessary for equation (53) 
to lead to meaningful results, the following restrictions are 
necessary : 

(1) The physical variable must be zero in the potential 
flow. 

(2) The physical variable must be homogeneous in the 
turbulent flow. 

(3) The physical variable (and its square) in the turbulent 
flow must be uncorrelated with the location of the front. 

(4) ddiere must be no mean value in the variable between 
turbulent and j)otential flows for the same value of //. 

The first condition immediately eliminates velocity fluc- 
tuations from this sort of treatment. This renders uncertain 
Townsend’s tuibulent energy application, mentioned above. 
However, vorticity fluctuation and turbulent shear certainly 
satisfy it, as may temj)erature or concentration fluctuations 
and heat or mass transfer, when these are ])resent. 

The second requirement is ])robably not satisfied by any 
variables in (lows with transport, including, of course, the 
commonest exam])le, shear flow. This follows from the 
fact that, evmn in spatial zones with 7 = 1 everywhere, there 
are gradients in all of the quantities which have been meas- 
ured. Consequently, the entire concept of Gt- as a function 
of ])osition in a shear flow must be semiquantitative at best. 

It seems unlikely that the third requirement is satisfied by 
all of the physical variables, but for most of them it may be 
close enough that equation (57) would be apjiroximately 
true. 
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Finally, the occuiTence of a mean value between potential 
and turbulent fluid must also de])end upon the particular 
ph^^sical variable under consideration. It certainly does 
occur for longitudinal velocity at the boundary of a turbulent 
shear flow. It certainl}" does not occur for z-component 
velocity in a shear flow which is two-dimensional in the 
mean with gradients all in the x- and 7/-directions. For 
many physical variables its occurrence or absence is not a 
priori obvious. In any case it can be handled by resorting 
to equation (57). 

P^or complex cases, when even equation (57) is believed to 
be inadequate, possibly because the variable is not zero in 
the potential zone, it is still possible to obtain statistical 
information on the signal structure within the turbulent 
bursts by laborious computational procedure for the oscillo- 
graphic trace. 

A more detailed question may be raised at this point as to 
the influence of intermittency upon the measured power 
spectrum of velocity fluctuation. Again the answer is 
certain to be simple if the four conditions listed above are 
satisfied. In that case, with probe signal a{t)=aT{t)b{t) , the 
autocorrelation functions of the three variables are related by 

^a{r)=\l/aj,{T)yph{r) (58) 

where r is time interval, dlie power spectra are simpi}' the 
P^ourier cosine transforms of the correlations and, since the 
transform of a product is equal to the convolution integral 
of the individual transforms, the three power spectra are 
relatecl bj^ 

-?’<.(«)= Fa^{nx)FJji — n^)dni (59) 


The first step is the strictly experimental process of ex- 
amining a new aspect of the degree of similarity to be found 
in the detailed structures of the various turbulent shear 
flows. 

The second has as its purpose the approximate verification 
of the rather crude hypotheses leading to prediction of the 
turbulence front behavior, that is, to equations such as (34) 
and (35). 

ROUGH-WALL BOUNDARY LAYER 


Pitting the experimental results on boundaiy-la.yer thick- 
ness by a simple power-law relation (see appendix), it turns 
out tliat, neglecting Rejuiolds number effects, 

bozb^Qceoz{x—Xof-^^^^-^ ( 60 ) 

numerically. 


5 ?«0.19(a:— in. 
5*^^0.13(:r— in. ^ 
(9«0.026(x-Jo)'' "' in.^ 


(60a) 


"Idle fitting of a power law to a set of points without origin 
involves two steps: (1) An origin must be chosen by trial and 
error to give the closest approximation to linearity on 
logarithmic graph paper, and (2) the “best’^ straight line 
must be drawn through the resulting plot. This procedure 
was also applied to the fitting of power-law approximations 
to the experimental data on <j{x) and Y{x). P^igure 21 
illustrates the degree to which a power-law fitting is success- 
ful. The latter quantities are then given by this '‘direct’’ 
measurement as 


(7«0.022(x-a:,)o*67±o.i 


( 61 ) 


where Fa{n) can be measured directly from the output of the 
hot-wire anemometer; F ^{n) = F c(n) , the spectrum of the 
Schmitt trigger output in the intermit tency-measuring 
circuit (see fig. 28) ; and Faj,{n) is the spectrum of a hypo- 
thetical homogeneous turbulence variable which should give 
the nature of the fluctuations within the bulges of the 
wrinkled front. 

Equation (59) is a PTedhohn integral equation of the first 
kind, readily solved in principle by P^ourier integral methods — 
which corresponds in effect to going back to equation (58). 

No attempt has been made to apply this relation because 
the experimental results ajipear too uncertain to merit such 
detailed manipulation. It is hoped, however, that such a 
study can be made in later shear-flow reseai'ch. 

APPLICATIONS TO PARTICULAR TURBULENT FLOWS 

Application of the foregoing general conce])ts and theoreti- 
cal predictions on the behavior of the turbulence front to 
particular turbulent flows involves two explicit aspects: 

(a) Comparison of directly measured a{x) and Y{x) with 
measured values of characteristic shear-layer thicknesses, for 
example, 0(x) ( oc5* oc5) in the boundaiy layer. 

(b) Comparison of (j{;x) and }"(-^), as computed from 
measured turbulence data, with directl}' measured values of 
(t{x) and Y{x). 


F«^0.14(x— in. (62) 

The "best” common origin is iTo= — 20 inches. 

Comparison of equations (61) and (62) with equation (60) 
shows that, witliin the precision of these experimental results, 
the turbulence front both progresses laterally and increases 
in amplitude at the same rate as the mean boundary-layer 
flow grows. The uncertainty range indicated is a crude 
estimate of standard deviation, not the maximum. 

Verification of equation (34) requires knowledge of both 
v'jF and (the transversal Lagrangian scale) as functions 
of X. On the basis of the v'j U measurements at corresponding 
])ositions across the boundary layer at four different j- 
stations (fig. 17), it is assumed for the sake of this calculation 
that v'ccUj, as dimensional reasoning and Laufer’s pipe 
measurements (ref. 15) also indicate. The Ur is obtained 
from the measurements of 6{x ) : Ur ocx~^‘^^ for very large 
values of x. 

Ihifortunately, there exist no measurements of Lagrangian 
scale in turbulent shear flows. However, the ratio of La- 
grangian to Eulerian scale has been measured as a 

function of v' Ljv for isotropic turbulence (ref. 18). These 
highly scattered measurements show L^IL to be a slowly 
decreasing function of v'Ljv. In order to estimate 
for substitution into equation (34) it is assumed that this 
variation holds roughly for shear flow. Further, there is 
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good ovideiu ‘0 that for a given shear Mow the Eulerian scale 
is ])ro])ortional to the eharaeteristie width of the shear zone; 
that is, in this ease, Lccb and the constant of ])ro])ortion- 
ality is taken from the smooth-wall boundary layer of 
Schubauer and Klehanotf (ref. 4) at a station where dl*ld.r =0. 
Their data give /> — (). 175. 

''.Vith this estimate of L/5, v' Ljv in the jiresent boundary 
layer goes from about to 500 in the primdpal test area: 
20 inches<j'< 1 1 0 inches. But over this range of v' Ljv, 
figure 54 of reference 18 indicates (by extra])olation) little 
change in Therefore, for ])urposes of the i)resent 

rough estimate, it is assumed that Li^ccLccd. With 
5a/" ®^ (ecp (()0)) and v' oz Vt the resulting 

theoretical ])rediction (eq. (54)) is 

a{x) oc (x—Xo)^'' (t)5) 


which agrees with the directly measured exponent (eq. (61)) 
])erhaps better than the accuracy of either measurement or 
t heoret ical a])proximat ion . 

Verification of equation (55) reciuires information only on 
^'(x) at corres])onding //-positions in the boundary layer. 
wSince this information is not yet directly available, one 


assumes the isotropic relation where v (x) has been 

measured and \{x) can be inferred by using the well-known 


isotroj)ic estimate^ 




and assuming 7. « 0.1 75 as 


before, losing the experimental value of 5(j) (eq. (60a)), 
this calculation gives 


r*(3’)oc (:r— Jo) 


(64) 


For com])arison with experiment this is next translated 
into y(x). Since V* is propagation velocity' relative to the 
fluid, one can write the ap])roximate relation 


dY^ y(V)+V* 

dx J](J) 


(65) 


which should hold for any reasonably fiat turbulence front. 
The term V is the mean fluid velocity in the //-direction. 
Equation (65) is approximate because (a) in some Mows the 
front is not very flat and (b) at y = Y the mean velocity of 
the turbulent fluid is_somewhat less than that of the non- 
turbulent fluid (i. e., U^). 

For the boundary layer, it is well known that 


cf5* 

dx 



( 66 ) 


which is easily shown from the definition of 5*(j). 


Since exj)eriments show that L\Y) ^ U{b) =V one can 
infer V{Y) =F(5), so that, for the boundary layer, equation 
(65) gives 


ilY^ _ V* r/5* 
dx ^ Ul'^dx 


(67) 


Since eciuation (60a) gives the exqierimental result 


dx 


0.08(j— j,)-^*39 


( 68 ) 


it is clear that the power-law a])])roximation to )^(j) will 
lie between 0.70 (if the F*-term dominates in eq. (67)) 
and 0.62 (if the 5*-term dominates). In fact, if the propor- 
tionality (‘onstant of e([uation (55) is determined from the 
data at a* =102 inches, 


F*^0.5Vrr 


(69) I 


If this is used with e(|uations (67) and (08) to predict 

f(j)«0.18(j-j.)'^-'^ (70) 

the agreement with the directly measured result, equation 
(62), is good. 

It should be remarked parenthetically that, although ade- 
quate measurements of t(j,//) are still not available on the 
smooth-wall turbulent boundary layer, an indirect verifica- 
tion of equations (54) and (55) follows from approximate 
agreement between the experimental (or anah^tically in- 
ferred) 5(j) and the ])redicted (t(x) and Y(x), using n'ason- 
ing like that ])resented in detail for the rough-wall case. 


TWO-DIMENSIONAL WAKE j 

The measununents of Townsend (ref. 10) in the plane wake j 
far behind a circular rod provide another case in which 
equations (54) and (55) can be checked against experiment. I 

From conservation of momentum and the assumption of , 
similarit}’', dimensional reasoning yields the experimentally j 
verified predictions that far behind the obstacle a turbulent | 
wake s])reads parabolicalh’- (5 oaxY^) and that the char- | 
acteristic mean velocity defect decreases ])arabolically 1 
{(r^ — I o) (ref. 50). This means that the plane j 

wake is a constant Reynolds numbei- shear flow and there- 
fore significantly simpler than, for example, the boundary 
layer. 

Far behind the wake-producing obstacle, where the fully 
developed wake is finally reached, the difi'erence between 
minimum and maximum velocity is so small that equa- 
tion (67) can be ajiproximated b}^ 


r/F^r* 

dx p 


( 71 ) 


wSince there are oidy two points in the full^" dev(‘loped j 
j^-range, it has not been ])ossible to determine em])irical 
])ower laws for 5, }", and a. The pertinent experimental j 
result is simply that within the experimental uncertainty 
the points in the fully developed range are consistent with | 
])arabolic growth for all three lengths. 

Verification of equation (54) again requires data on ^ (j) 

L> 

and Lf^{x). As can be antici])ated for a constant Reynolds 
number flow, the root -mean-square turbulent velocities are , 
])ro])ortional to the characteristic mean velocity (mean 
velocity difference in the wake) so that r'ocj"^''^. Tn this 
asymptotic state, the mean velocit.y differences are all 


® The constant of proportionality is obtained empirically from reference 29. 
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small — so that Constant Further- 

more, constant Reynolds number implies Lagrangian scale 
proportional to Eulerian scale (/>lW ozL{x))^ and with the 
general assumption of Lccb it follows that P]qua- 

tion (:H) then gives as prediction for the variation in wrinkle 
amplitude of the turbulence front 

(j{x) oc.r^^^ (72) 


For the evaluation it is again assumed that J'ocy'/k- 

With oz{x—Xo)~^ and \ocL ca\i 2 Cc{x—Xo) equation (35) 
gives 

oz{x-xo)-^ (79) 

and the comparison with experiment can be made by using 
equation (79) to predict R{x), merely replacing Y by R in 
equation (65). Instead of attempting a detailed proper cal- 
culation only a rough estimate was made by assuming 


in agreement with the directly measured result, in the simi- 
larity (large x) zone. 


For the comparison of equation ( 35 ) 

\ with experiment, no 

data on are available. As in the boundary layer, it will 

therefore be assumed that ocr'/X. 
equation (35) predicts 

Since v' ccx and 


(73) 

whence 


YozxY^ 

(74) 


again in agreement with the directly measured result. 


ROUND JET 

Since fairly detailed turbulence data were already avail- 
able for the case of the round turbulent jet (refs. 2 and 31) 
entering fluid at rest, intermittency surveys 7 (j',r) have been 
made during the course of this investigation to provide 
further experimental check on the proposed physical picture 
of the turbulence front. 

From conservation of momentum and the assumption of 
similarity dimensional reasoning yields the experimentally 
verified predictions that far from its source the round turbu- 
lent jet spreads linearly {r\i 20 cx) and that the chaz*actei-istic 
mean velocity decreases hyperbolically (f oOcj^“^) (ref. 30). 
Thus, the round jet is another constant Reynolds number 
flow and therefore relatively simple. 

The new measurements made in the course of this study 
(fig. 24) give as power-law approximations with the “best’’ 
common origin. 



(rf ‘ 

\ I.00±0.05 

’) 

(75) 

77oc^ 

1-^ 

^0.88±0.05 

(70) 

aoc( 


V 1.06±0.05 

(77) 


which may all be taken as linear within the experimental 
uncertainty. 

Frevious measurements have shown y'/f to be constant 
and independent of x at corresponding radical positions in 
the jet. Furthermore, the constancy of RcAmolds number 
again permits the inference that LlocL. With the assump- 
tion that LoeXi/Of equation (34) predicts 


■OC{x—Xo) 


in reasonable agreement with equation (77). 


V(x,R) ocU{XyR) (xU(x,0) 
Then the prediction is 


R(x) oc(x-Xo) (80) 

in reasonable agreement with equation (76). In fact, a be- 
lief in full similarity for constant Reynolds number shear 
flows suggests that equation (80) may be more nearly correct 
than equation (76). 

INTERMITTENCY AND MEASURED MEAN QUANTITIES 

As pointed out earlier in the section on ‘Tnfereiice of 
Turbulence Properties From Intermittent Signal,^’ there 
seems to be only a restricted likelihood of extracting from 
the measured statistical characteristics of the intermittent 
signal respectable quantitative results on the statistical 
properties of the turbulent flow in the convex bulges of the 
turbulent front. Probably the broadest obstacle to simple 
physical interpretation of results computed from eci nation 
(53) or (57) is the lack of homogeneity within a fully turbu- 
lent zone supporting transfer. 

Nevertheless, it seems worth while to present, for some 
fluctuating variables which are zero in the potential fluid 
(i. e., satisfy the first requirement), the results of applying 
these two operators. 

No detailed quantitative information has yet been obtained 
on the mean- value jump for any physical variable. A rough 
check from u{t) oscillograms in the rough-wall boundary 
layer where 7 — 0.4 indicated that the jump in longitudinal 
velocity was about 5 to 10 percent of f „• 

Unfortunately, this still does not permit calculation of uY 
because all velocity fluctuations violate the first condition; 
that is, they are nonzero in the potential flow. 

For UTorticity, which does satisfy this first condition, no 
mean-value jump is observable on the oscillograms. This is 
not surprising since this j-component has no corresponding 
mean vorticity in this flow field. It may be anticipated 
that the ^-component f will be found to have a jump, if and 
when it is measured. A plot of = is given in figure 30. 
To insure ?/-coordinate consistency, this particular 7 ( 7 /) has 
been measured with the vorticity meter as sensing element. 
It does not differ appreciably from 7 ( 7 /) as determined from 
the differentiated signal of a i/-meter. 

In the outer part of the intermittent zone p /7 turns out 
to be roughly constant, leading to the possible conclusion 
that is relatively homogeneous in the turbulent fluid. 

The Reynolds shear stress —puv has been inferred to be 
zero in the potential field outside a turbulence front. There- 
fore, it may also be interesting to estimate —pUrVr. 


(78) 
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Fk: RE 30. — Di.stribution of moan-sciuan' vorticity fluctuation divided 
by intennittcMicy in boundary layer at x=102 inche.^. 

Using a ropreseiitalion like equation (54) for u{i) and v{t) 
separately and assuming (a) no m ean jump in v{t), (h) 
={), and (c) h'^UTi-T it follows that 

(81) 

Townsend (ref. 28) has plotted uvjy^^ versus y without 
att('mpting a justification. 

F'igure 31 shows the residt of applying equation (81) to 
the measured Reynolds shear stress in the rough-wall 
boundary layer and in Townsend’s plane wake (ref. 10). 
TIh' nonconstancy of uv/y can probably be attributed largely 
to nonhomogeneity within the turbulent field. 

In concluding this section it may ho remarked that, if an 
existing nonhomogeneity for any variable in the turbulent 
part of the field depends only upon distance in from the front, 
a first-order estimate of its effect can be made by computing 
the average value generated at a fixed point by random 
motion of a “fixed” pattern like that in sketch (e). 



STATISTICAL DESCRIPTION OF TURBULENCE FRONT 

The position of the turbulence front I is a random 
variable stationary in time and nonstationary in x. The 
pur])ose of this section is to report some further measure- 
ments which have been made on its statistical properties, 
especially those of Y(f) for a fixed value of x. Earlier 




(a) Boundary layer at x=102 inches. 

(b) Plane wake at a: — 800 inches. (Data from nT. 10.) 

Fh:ure 31. — Disi ributions of lUwnolds shear stress divided by inter- 
mit tency. 

sections have emphasized its statistical variation with Xj 
particularly through r(-^) ^Did the standard deviation c(x). 

It is of course possible for Y to be a multiple-valued func- 
tion (see, e. g.. fig. 1), but in most (lows the occurrence of 
multiple values apj)ears to be sufliciently rare that a discus- 
sion predicated upon a single-valued Y is applicable with 
good accuracy. This is especially true for the boundary 
layer, where turbulence levels tend to be appreciably lower 
than, for exani|)le, in jets entering a still medium. This 
conceptual restriction to single-valued Y, exercised through- 
out the report, will be justified empirically for the boundary 
layer by showing that the average wave length is considerably 
greater than the average wrinkle amplitude. 

As a stationary random function )\t) is susceptible of 
quantitative statistical description in various ways, not all 
indei)cndent. Perhaj)s the two most common mutually 
independent functional representations for such variables 
are the autocorrelation function (or its F^ourier transform, 
the power spectrum) and the probability density (or its 
F\)urier transform, the characteristic function). Usually 
the lower order moments of the density and spectral func- 
tions, which have simple i)hysical interpretations, are the 
most easily measured statistical properties. 
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The problem of acquiring detailed statistical information 
on Y(t) is novel in the sense that nowhere in the experiment 
is there a signal which is simply proportional to the stationary 
variable under study. Therefore, the conventional statistical 
functions (above) are not readily measurable by standard 
techniques. It is fortuitous that the dissimilar character of 
the fields on opposite sides of Y(t) gives such a convenient 
method of measuring probability density. However, the 
autocorrelation or power spectrum apparently cannot be 
directly measured, and tlierefore other direct statistical data 
have been sought, in particular, the probability density of 
pulse lengths, actually the statistical measure of the times 
between successive occurences of any particular value of the 
])rimary variable Y(t). 

A challenging problem in the theory of stochastic processes 
is that of relating (if possible) these densities to the more con- 
ventional statistical measures. Up to the present time, only 
a few fringe results seem to have been obtained by workers 
in the field; these will be mentioned in approjiriate context. 

PROBABILITY DENSITY OF Y{t) 

As has been pointed out in the ^Hntroduction’’ (eq. (D), 
the intermittency factor y{y) is simply the distribution 
function of Y{f) and, therefore, dy/dy is its probability 
density. 

Calculation of dy/dy shows that, except in the two tails of 
the function, it is remarkably symmetrical. Furthermore, 
the physical picture given here of front wrinkling as pri- 
marily a (Lagrangian) turbulent diffusion phenomenon then 
suggests a check to see how nearly dy/dy approximates a 
Gaussian function, since studies of scalar diflhision in iso- 
tropic turbulence have shown a closely Gaussian density. 
Figure 32 shows this check. It includes typical plots on 



Gaussian paper of y{y) for the boundary layer and for Town- 
send’s wake, as well as y(r) for the round jet. Clearly all 
three distributions are Gaussian within the experimental 
precision except in the tail regions. 

Deviations from s^unmetry must, of course, occur at tlie 
tails since the boundary conditions on the two sides are 
vastly different. 

wSince the nearly Gaussian character of dispersion in iso- 
tropic turbulence is still unexplained theoretically, it is not 
to be expected that this much more complex plienomenon 
can be clarified at present. Also, it must be emphasized that 
even in the former case it is not necessarily true that the 
probability densities are precisely Gaussian; the current con- 
clusion is only that a Gaussian curve fits the data as closely 
as present experimental techniques produce data. Very 
likely it is the deviations (however small) which, when 
measured, will shed more light upon the central property of 
turbulence, the nonlinearity. 

Batchelor (ref. 32) has pointed out that the Gaussian dis- 
])ersion ])attern observed at very large distances downstream 
from a contaminant source in a turbulent flow may be simply 
a consequence of the central limit theorem,^ since the relative 
position of a fluid particle a long time after tagging may be 
regarded as the sum (time integral) of a large number of 
small dis])lacements, which are at least uncorrelated for mod- 
erate intervals if not exactly statistically independent. In 
fact, if this reasoning does apply, it is doubl}^ effective: 
l^article displacement, the jirincipal variable, is itself the 
integral of particle velocity, so that the long-time displace- 
ment is the sum of a collection of sums. 

Apparently, the central limit theorem has not been ex- 
tended to integrals of continuous random variables, but some 
pertinent work has been done by Kac and Siegert (ref. 33), 
who sliowed mathematically that passage of a particular 
skew (probability density) random signal through a low-pass 
filter reduces the skewness. This prediction has been experi- 
mentally verified by Jastram (ref. 34) and by Iribe (ref. 35) 
A low-pass filter is, of course, qualitatively equivalent to 
integration. 

PROBABILITY DENvSITY OF PULSE LENGTHS 

Experimental results. — Fi-om a sketch of Y{t) as a sta- 
tionary random variable, it is easily seen that the intermit- 
tent signal from a fixed probe provides a direct means of 
measuring the statistical distribution of the time intervals 
between successive occurrences of any particular value of F. 



Figure 32. — Typical intermittency di.stributions for boundary layer, 
jet, and wake plotted on Gaussian probability scale. 


7 In effect, this states that the sum of a number of statistically independent random variables 
approaches Gaussian character as the number increases without limit (provided that no 
finite group dominates the sum). 
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From the fluid-moclianical point of view this gives a eoiiven- 
ieiit measure of the wave lengtlis of the mountains and valleys 
in the turbulenee front. The laminar superlayer is thin 
enough to be considered a discontinuity in all of this analysis. 

If 1\ is the duration of the probe in a turbulent zone and 
7^2, the duration in a nonturbulent zone, figure 29 gives the 
probability densities p\{Ti) and jhiTo) at three different 
values of transversal position, that is, tliree different values 
of the intermittency factor. 

By definition (of probability density), the curves in figure 
29 are normalized to unit area. A check on their accuracy 
is given by tlie more or less obvious condition 


win ‘re 


JI\_ 

7 \+ 7 ; 


7 




ipi(7\)(/7\ and 



The terms 1\ and are average j)ulse durations in units 
of time and are functions of y or, alternatively, of 7 since 7(y) 
is monotonic. 

The computations from figure 29 give (for 5 — 3.5 inches): 


y/5 

y 

directly 

measured 

Ti, sec 

To, sec 

7 _ 
from T\ 

and To 

li, in. 

h, in. 

0. 72 

0.75 

0.0155 

0. 0000 

0. 72 

7. 0 

2.7 

0. 85 

0.50 

0.010(3 

0. 0082 

0. 50 

4.8 

3.7 

0.98 

0. 25 

0. 0009 

0.0132 

0.34 

3. 1 

5.9 


where U = T^1\ and h = l^o.T 2 are approximate measures of 
the sj)atial extension of the average intervals in this x vicinity, 
ddiis interpretation of the /’s as average intercept lengths for 
the random variable Y{x) gets increasingly accurate as the 
velocity fluctuation level decreases. This time-space trans- 
formation is, in fact, identical with that first proposed by 
Taylor for an isotropic turbulence (ref. 36) and discussed in 
moi*e detail by others (refs. 37 and 18). 

A comparison between^ (/i +/2) for 7=0.50 and the stan- 
dard deviation a of }^i(0 at the same x-station gives a rough 
measure of the flatness of the wrinkled turbulence front. 
For this particular station in the boundary layer, 


/1 + /2 


«0.13 


(83) 


which indicates a rather flat front, as assumed in the eai*lier 
theoretical discussion on the j)ro])agation of tlie laminar 
superlayer. 

Inspection of figure 29 shows the following traits of the 
data ; 

(a) The points are rather scattered. 

(])) For 7 -0.50, Pi and yn show an aj)preciable difference. 

(c) The 7=0.25 and 7=0.75 cases, which might be ex- 
pected to have identical curves with reversed labels, show 
this character qualitatively, though not accurately. 

I’roperties (b) and (c) can ai)|)arently be attributed chiefly 
to the shortness of oscillographic samples;^ therefore, the 


curves in figure 29 have been labeled with the 7’s actually 
given by these short samjdes, and the apparent discrepancies 
(b) and (c) are qualitatively explained. In other woi'ds, a 
short sample with actual 7=71 drawn from an infinite record 
with 7=72 can be expected to show other statistical proper- 
ties resembling those of an infinite record with 7=71. 

Two other sources of uncertainty in the data of figure 29 
are (1) the natural uncertainty of measurement in the 
presence of noise, even with perfect equipment, and (2) im- 
perfections in measuring equipment and technic|ues. 

The first of these difficulties affects all intermittency meas- 
urements and is basically insurmountable. Of cours(*, the 
noise level could be reduced somewhat and, under simplifying 
statistical assumptions on both noise and signal, some esti- 
mate of the effect could be made. 

The second difficulty probably affects pi and /n measure- 
ments more seriously than direct 7 measurements. For ex- 
ample, su|)pose that the measuring process misses a sizable 
number of the shortest turbulent bursts. This fault will 
scarcely affect the directly measured 7 since these contain 
only a small part of the total number of pulses to be co.inted 
(except for y<^A). On the other hand, this fault will not 
only change the character of P\{7\) for small values of Ti but 
also will change the level of yn(T2) for large values of T 2 , since 
the very short turbulent bursts subdivide long potential 
bursts into shorter ones. Hence, this fault will serioudy 
affect and, therefore, 7 as computed from 7\ and T 2 ' 
Precisely this fault is observable on the oscillographic traces. 

Other defects similarly observed are the (relatively in- 
frequent) missing of short potential bursts and the occasional 
overhang of the trigger signal beyond the duration of a tur- 
bulent burst. The last of these faults affects the direct 7 
measurement as well. 

An obvious way around some of these difficulties is the 

direct use of ^{t) cr^ {{) oscillograms to compute yq and y>2. 

To some extent this was done, and the extreme tediousness 
of tliis method is exactly why the samples processed are so 
short. 

This inadequate sample length (fault (3)) most seriously 
affects the results in the large and T 2 ranges. The relative 
seriousness of this limitation for long versus short pulses is 
not given (as might be guessed at first blush) by the ratio of 
sample length to pulse length but by the ratio of sample 
length to the inverse of the frequency of occurrence of the 
])articular length of pulse (actually a small range) in question. 
For example, in a 3-second oscillographic sample, the i)oints 
on the tails of p\ and pi may repi’esent as few as one or two 
actual occurrences. Witli this in mind it can be concluded 
that tlie agreement between values of 7 obtained via Ti and 
72 and values of 7 directly measured is surprisingly good. 

It would be interesting to know whether pi and pi approxi- 
mate exponential distributions for large values of Ti and 7V 
However, the uncertainty of the j)oints in just tins range is 
so great as to render such a (juantitative question unanswer- 
able. Some v(‘iy indirect evidence via the power spectrum 
of the Schmitt trigger output for 7=0.50 will be discussed 
in a following section. 

Since the small Ti and T 2 ranges of pi and pi are quite un- 
certain (i.e., for bursts shorter than 2 milliseconds), some 


* About 3 socoiuls, as coniparocl with the 2 niiniitos used in obtaining the 7 *s directly. 
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qualitative analytical consideration of the anticipated be- 
havior in this range is in order. These short segments arise 
whenever the hot-wire ])asses just below a local maximum 
(for turbulent fluid) or just above a local minimum (for po- 
tential fluid) in 1"(0- 

The variable V(f) must he differentiable (since it occurs in 
a continuum); therefore, its extremes have horizontal 
tangents. Thus, a Taylor series expansion of )"(/) about any 
local extreme t = starts witli a term proportional to {t — tmY. 
The limiting behavior of pi and ]) can thus be obtained by 
considering a parabola C=I)~ as in sketch (g). Tlie problem 
is then as follows: Suppose C has a flat probability density 
Vc{C)Y what is the probability density po of 


follow easily from Rice’s result. Of course, in the particular 
case when successive intervals are statistically independent, 
the occurrence numbers have a Poisson density, and the 
interval lengths have a sim])le exponential probability 
density. 

A more directly a])j)licable result, apparently due to Rice 
(ref. 38), relates the expected rate of occurrence of any 
particular value of a Gaussian variable I{t) to the probability 
density of the variable and the autocorrelation function 
behavior in the vicinity of zero: 



In general, if C—C{I)) is unique, 

VdU>)=Vc{C) (84) 

whence, for this particular problem 

p^{l)) ozD (85) 

It follows from this calculation that the probability density 
of the intervals between successive occurrences oj any jmrticular 
value of a dijferentiable random variable must start out {a) 
from the origin and (6) linearly. Specifically, P\{Ti) and 
7 ^ 2 ( 72 ) must behave in this fashion, even though the measured 
curves do not all show this tendency in the range covered. 

By reasoning similar to the above it is obvious that for a 
continuous but nondifferentiable variable (corresponding to 
pointed but uncusped extremes) the corresponding density 
starts out at a finite value. 

Status of random-variable theory. — ^The mathematical 
problem of relating the probability density of the intervals 
between successive occurrences of any particular value of a 
continuous random variable to the ordinarily more accessible 
statistical functions (probability density of the primary 
variable, power spectrum, etc.) has apparently not been 
solved, even for a Gaussian variable. 

Rice (ref. 38) has deduced the probability of a zero of a 
Gaussian variable 7(f) in an interval {U-Yto) XhY-to + dt) 
when there is a zero at U. However, the probability density 
of intervals between successive zeros (or successive occur- 
rences of any other particular value) does not appear to 

» The very small range to be studied, that Ls, just the immediate vicinity of an extreme, 
permits approximating any small segment of a finite probability density by a constant value. 


where \f{T) is the nonnormalized autocorrelation function 
I(t)I{t-\-r) and a prime indicates differentiation. The 
proof of equation (86) requires also that 7(0 and 7'(f) be 
uncorrelated — which is automatically satisfied for a station- 
ary variable. 

However, it must be emphasized that the pristine simplicity 
of this theorem is dependent upon the restriction to a 
Gaussian variable. Two of the seemingly inexhaustible 
number of fortuitous properties of the Gaussian probability 
densit}^ are: 

(a) If a variable is Gaussian, so is its derivative. 

(b) If two Gaussian variables are uncorrelated, it follows 
that they are statistically independent. 

Without tliese built-in conveniences, it seems likely that 
sucli a theorem could be deduced only with the general 
assumptions tliat the variable and its derivative are statisti- 
cally independent. 

For the expected rate of zeros, equation (86) reduces to 
(ref. 38) 



Equation (87) has been used by Liepmann, Laiifer, and 
Liepmann (ref. 39) to measure the microscale in a decaying 
isotropic turbulence. It will be used here to obtain yp" 
for the turbulence front Y{t). 

Measurements have also been made of the average rate of 
occurrence of the values of Y{t) corresponding to 7=0.25 
and 0.75 in the intermittent zone of the rough-wall boundary 
layer. In figure 33, the three experimental points are 

y-y 


compared with equation (86). The 


location of the 


three experimental points has been chosen according to the 
value of 7 of the short samples (from which the iV's were 
measured; see the preceding table) rather than the true 
physical locations of the probe. The agreement is better 
than can be expected with the uncertainly of the measure- 
ments and therefore fortuitous. The number given for the 
rate of occurrence of zeros, A^o~108 per second, is inter- 
polated along the Gaussian curve. 

Measurements of the probability density of zeros in the 
fluctuating part of (he signal from a Jiuinan voice have been 
reported by Davenport (ref. 40). 
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where Yi = Y — Y. Obviously \l/{0) = a~ and, since A^o ^ 
are the measured quantities, equation (87) is written 


where 


,/."(()) = -7tVAV 


2 

A^0 = ^ I Tp zi'i'os/sec 


(90) 


(90a) 


For the rough-wall turbulent boundary layer at x ==102 
inches, o- = 0.55 inch and A^c=108 zeros per second, so 


Fir. r RE — Frociuency of occurn'iice of zero and two other particular 

values of Vi(0 in lioundary layer at x=102 inches. Solid curve is 
that of a strictly CJaussian variable. .Vo=108. 


POWER SPECTRUM OF SCHMITT TRIGGER OUTPUT 


As indicated in figure 0 , the output of the Schmitt trigger 
is in principle a random flat-top signal which is on whenever 
the ])robe is in turbulent fluid and oft* whenever it is in poten- 
tial flow. Obviously the statistical properties of this signal 
must have some relation to those of the primary variable 
Y{t)y and therefore two convenient properties have been 
measured. The first is the probability density of pulse 
lengths, tops and bottoms separately; these are, of course, 
just pi(Ti) and 7 ^ 2 (T 2 ) (fig. 29). The second is the power 
s])ectrum of the trigger output, measured at the same loca- 
tions as the densities (fig. 28). 

Tlie three spectra have the same general shape, with power- 
law decrease for high frequency as indicated in the figure. 

It might be expected that a relation should exist between 
the ])ulse-length densities of any flat-top signal and its power 
s])ectrum, but a search of the literature has uncovered no 
such analytical results except in special cases, one of which 
is used below. 

The simplest of the three signals is that corresponding to 
-y^ = 0.50, and in figure 28 this power spectrum is seen to agree 
closely with that for a ^'Poisson type'’ flat-top signal (see, 
e. g., ref. 38) : 


F(n)oz 


M 

A/-+7t2/P 


(88) 


whei’e M is the average number of jumps per second and n is 
cyclic frequency. For this application and 7 = 0.50, 
the average number of zeros per second in Y — since the 

distribution of zeros cannot be truly Poisson. 

The very good agreement in figure 28 implies only that in 
this case P\(Ti) and poiTo) could be exponential away from 
the origin, even though the directly measured data are too 
uncertain to permit any estimates. However, no assertion 
can be made, since tbe —2 power spectral decrease is charac- 
teristic of most signals with ‘hliscontinuities.” 


■p" (0) = 3.4X10^ sq in./sec- (91) 

A corresponding characteristic length mathematically equiva- 
lent to the dissi|)ative scale (microscale) in turbulence can 
be deduced by the time-space transformation mentioned 
earlier: 



For this particular case. 


Xy^«1.9in. (93) 

which is a bit smaller tban li and L in the preceding lable 
for 7=0.50. 

For low turbulence levels, one miglit expect the ([uantity 

= Xy' to be of tlie order of the Lagrangian spatial microscale 

U 

(ref. 18), which is roughly ecpial to the Eulerian 
microscale X over a wide range of I\\ in isotropic turbulence 

(ref. 18). In this case, Xy^^j=:^0.09 inch. Tliis is the same 

order as X in the neighboring turbulence. Since A\ for this 
turbulence is roughly 70, which (in isotropic turbulence) 
gives X^«1.5X, the conclusion here is that 

= X,-«y'X, (94) 

or, since the Lagrangian time mocroscale of the 

neighboring turbulence is given roughly by 

Equation (90) gives only the vertex curvature of the auto- 
correlation function. ILnauise of the Gaussian character of 
Y{t), it is possible to estimate the entire Pir) from the spec- 
trum of the trigger output. It has been shown by North 
(see ref. 41) that the autocorrelation function of a strongly 
clipped Gaussian variable is simply related to the autocor- 
relation function of the variable itself: 


AUTOCORREL.\TION FUNCTION OF ¥(() 

The approximately Gaussian character of Y^ Y permits 
application of equation (87) relating the zero occurrence rate 
and tlie autocorrelation. For this jiurpose the nonnormal- 
ized autocorrelation is defined by 


P^(r) 
Pc (0) 



- 1 



(96) 


A strongly clipped variable is just a flat-top signal which 
changes sign whenever the primary variable passes through 
zero — which exactly describes the relation between the trigger 
output and the primary variable Y(t). 


P(T) = Yi(f)Yi(t + r) 


(89) 


FREE-STREAM BOUNDARIES OF TURBULENT FLOWS 
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Since, as shown by Wiener (ref. 42), the autocorrelation 
function of a stationary random variable is just the Fourier 
cosine transform of its power spectrum (and vice versa). 


^c(r)=J^ Fein) 

00 

Fc (n)=4 J ype { t ) i^os 2Trn 


cos 2Trn T (In 
rdr 


(97) 


the autocorrelation of the trigger output is computed from 
the measured power spectrum. The good agreement of 
Fc(n) with the form in equation (88) permits using a simple 
exponential for ypei^) (I'f'f- 38): 




Then, equation (94) gives the autocorrelation function of 
the turbulence front location: 


which is plotted in figure 34. The vertex osculating paral)ola 
corresponding to yp"{0) as given by equation (87) and the 
directly measured zero occurrence rate are drawn in for com- 
parison. The former parabola should give the origin be- 
havior of ypir) mmre accurately than equation (99). 



Figure 34. — Autocorrelation function of turbulence front location 
as a function of time in boundary layer at x = 102 inches. 


As should be expected, the calculation of i/'"(0) for equa- 
tion (99) gives 


^"(0) 

m 


= -ir\P 


( 100 ) 


identical with equation (87), if d/-^A^o- 

In fact, it found experimentally that \I^2Nq for this in- 
vestigation. This is not surprising since the differentiability 
of Yi{t)y whose zeros give the square-wave jumps, leads to 
a considerable deficit of short pulses as compared with a 
truly Poisson sc[uare wave (see the section ^‘Experimental 
results'' under “Probability Density of Pulse Lengths"). 

Of course, the power spectrum of Yi{t) could be calculated 
by taking the Fourier cosine transform of \p{r) but the data 


are sufficiently inaccurate that further manipulation scarcely 
seems worth while. 

Other characteristic lengths of the wrinkled turbulence 
front can be estimated from the integral of ypir)^ mathemati- 
cally analogous to the integral scale of turbulence, but these 
may be less pertinent than, for example, U and U, the average 
pulse lengths: 



whence, 

i(0)Jo 

which turns out to be the same order as /i and I 2 . 
Alternatively, 


botli values being for the rough-wall boundary layer at 
j:= 102 inches. 


CONCLUDING DISCUSSION 

From the analytical and experimental results reported 
here on the problem of the relatively sharp instantaneous 
front separating turbulent fluid from nonturbulent fluid (as 
at a free-stream boundary), the following new conclusions 
are drawn: 

1. The nonturbulent region is a field of irrotational fluc- 
tuations. 

2. The front separating turbulent from potential flow is 
actually a very thin fluid layer in which viscous forces are of 
primary importance. The role of this “laminar superlayer" is 
the propagation of vorticity (both mean fluctuating) into 
the potential field. It is m.aintained thin by propagation 
relative to the fluid and by the random stretching of vortex 
lines in its local vorticity gradient. 

3. The common occurrence of contiguous rotational and 
irrotational velocity fluctuation fields underscores the useful- 
ness of confining the word “turbulent" to random rotational 
fields only. 

4. The rate of increase of wrinkle amplitude of the turbu- 
lence front can be roughly predicted in terms of a Lagrangian 
diffusion analysis, using the statistical properties of the 
turbulence in the fully turbulent zone. The actual estimate 
is given by equation (34). 

5. By dimensional reasoning and, independently, through 
a model of the laminar superlayer, the thickness of the super- 
layer can be estimated. The simplest approximation is 
equation (51), giving a thickness of the same order as the 
Kolmogoroff (minimum) turbulence length. 

6. The propagation velocity F* of the turbulence front is 
taken by dimensional reasoning to be proportional to V^- 
This is roughly verified by experiment. 

7. The downstream rate of growth of the turbulence front, 
as measured by standard deviation (j{x) and transversal 
position Y{x)j is found to be proportional to the shear-zone 
thickness, within the experimental precision, for plane wake, 
round jet, and rough-wall boundary layer. This is shown 
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iiuloiieiulcntly by direct experiment and by application of 
the results outlined in the i)revious paragrai)lis. 

8. The probability density of the turbulence front location 
at any fixed downstream station is (iaussian within the 
precision of the measurements everywhere except at the tails. 
This is found experimentally for all three types of turbulent 
shear flow studied. 

9. The probability density of the pulse lengths in the 
intermittent signal deviates strongly from the simple ex{) 0 “ 
nential type, presumabl}^ because V(t) is differentiable. 

10. The autocorrelation function of Yiit) for the boundary 
layer is found very indirectly from experiment to be as 
shown in figure 34. 

It seems likely that the ])reseiice of the turbulence front 
with its attendant detailed statistical properties will have 
to be included in basic research on turbulent shear flows 
with free-stream boundaries. It is not quite so clear that 
it must be explicitly included in semiemjiirical engineering 
estimates concerned only with overall transfer; so far no 
case has been encountered in which the front grows at a 
rate distinctlv different from the gross shear-layer growth. 

It appears that at present this new physical picture 
introduces at least as many new questions as it gives expla- 
nations of older observations. Insofar as it is concerned 
with a boundary condition, it tells nothing about transport 

phenomena within a turbulent region. Yet, since the 

wrinkle amplitude a{x) and transversal travel V{x) of the 
turbulence front appear to be governed by (or related to) 
[iroperties of the contiguous turbulence, any gross assump- 
tion on these variables implies consequent relations among 
the turbulence properties. 

It should especially be pointed out that the present 
investigation does not appear to shed any light on the 


characteristic difference between transport rates of vector 
(momentum) and scalar (heat, mass) properties. In fact, 
since it is concluded that no mean momentum can be trans- 
ported beyond the turbulence front it appears that (for 
laminar Prandtl and Schmidt numbers not very much smaller 
than unity) the front should apply equally well to heat or 
chemical composition. Oscillographic observations (not 
mentioned in tlie body of the report) in a hot jet show a 
temperature fluctuation intermittency, presumably coinci- 
dent with the vorticity intermittency. If this infereiu'e is 
true, then the vector versus scalar transport rate difference 
will have to be explained in terms of jiroperties of the en- 
tirely turbulent region. 

Interesting speculations in this direction have been made 
by Townsend (ref. 10), who suggests that momentum is 
largely trans})orted by relatively high wave number fluctua- 
tions while heat is transported by l)oth low and high wave 
number fluctuations, that is, by jet convection and by 
gradient diffusion, respectivel}L However, there are two 
dubious minor postulates in his analysis (mentioned here 
in the section ‘‘Inference of Turbulence Properties Prom 
Intermittent Signar' and at the end of the section “Laminar 
Superlayer’’) and also he has not clarified the jirincipal 
assumption vis-a-vis the known fact that the shear correla- 
tion uv appears to get ever increasing contributions toward 
the low wave numbers (ref. 43). P^inally, his inference that 
the lateral jets (bulges) convect little longitudinal momentum 
appears to be in contradiction to the fact that the inter- 
mittent velocity signal shows an a|)preciably lower mean 
in the turbulent segments tliat in the potential ones, as 
seen in figures 9 and 15. 

The Johns Hopkins University, 

Baltimore, Md., January 20, 1053. 


APPENDIX 

GROWTH OF ROUGH-WALL BOUNDARY LAYER 


Although the growth of tiirl)ulGiit boundary layers with 
zero static-pressure gradient is better approximated by a 
logarithmic function (ref. 11), the exploratory purposes of 
this investigation are satisfied by the simpler and less 
accurate power-law treatment. 

The momentum integral relation for tiirl)ulent ])oundary 
layer with zero static-pressure gradient can be written 
approximately as (ref. 11) 




(Al) 


The following rough assumptions are made: 

(a) Simple geometrical similarity in mean velocity profiles: 




(b) “Fully rough” wall conditions: 
Urh 

Therefore, 


>100 


r„<xp[^7(/i)]^ 

where h is effective roughness height, 

(c) Power-law velocity profile: 


■KIT 




From assuin])tions (a) and (c) 

r(h)_/ii 




whence the second assumption gives 




(A2) 

(A3) 


Since doc8, substitution of ecpiation (A3) into equation (Al) 
gives 


ax 

(A4) 

for h = Const ant . Therefore, 


1 

(Ao) 


p](luation (Ao), a simple power law, permits approximation 
to the actual boundary-layer growth with accuracy adequate 
for the present investigation. 

In fact, since both m and the exponent in equation (A3) 
have been measured independently, there is opportunity for 


an experimental check on the accuracy of the present crude 
ai)])roach: Mean velocity profiles (fig. 11) give 1/3.5. 
Tlierefore, the analysis predicts 

6cc(x—Xoy-^^ 

whereas measurements of boundary-layer growth (fig. 12) 
give 

It should be pointed out that boundary layers in general 
cannot have simple geometrical similarity because their 
characteristic Jieynolds numbers increase with x. 

This particular “rough-wall” boundary layer is fully rough 
all the way downstream (from x=0 to x=\{)2, falls 

from 200 to 145), if the peak-to-peak height of the corruga- 
tion is interpreted as h. 
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